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Resumen: En esta primera parte, después de resumir y reformular la teoría de bidominio 
de la electrocardiología, se emplea un modelo matemático simplificado de la excitación 
del miocardio a través de electrodos de marcapaso para estudiar, justo hasta el umbral, la 
dinámica del sistema electrodo-miocardio. Se introduce una región de influencia del 
electrodo de marcapaso sobre el miocardio, para poder plantear un problema mixto, ho- 
mogéneo de Dirichlet y no homogéneo de Neumann a partir de las ecuaciones no lineales 
del modelo de bidominio, de modo de describir el efecto polarizante del electrodo. Un 
enfoque basado en el análisis modal no lineal se emplea para transformar las ecuaciones 
de campo en un sistema de ecuaciones diferenciales ordinarias no lineales en las amplitu- 
des de modo. El valor umbral de la perturbación que causa inestabilidad se estima em- 
pleando una versión modificada del método de Eckhaus. Se identifican patrones espacia- 
les de polarización de las membranas que corresponden a una familia de estados umbral 
en el espacio de estados construido a partir de las amplitudes de modo. Se obtienen las 
formas y los tamaños de las regiones liminales (que deben ser despolarizadas por enci- 
ma del umbral uniforme de membrana para que pueda generarse un potencial de acción 
propagado). Se deducen fórmulas para curvas intensidad-duración (I-D), obtenidas me- 
diante pulsos catódicos de corriente controlada, incluyendo fórmulas analíticas para las 
constantes de tiempo (cronaxias) y reobases, en términos de parámetros geométricos y 
electroquímicos del electrodo, de la distancia entre el electrodo y el miocardio ventricular, 
de las propiedades eléctricas del conductor de volumen de los tejidos, y de la distribución 
local de los parámetros electrotónicos y de excitabilidad del músculo cardíaco. Se expli- 
can varios resultados experimentales (las regiones liminales y el comportamiento de los 
parámetros de las curvas I-D) y se predicen algunos fenómenos nuevos (relacionados con 
las regiones liminales, las familias de estados umbral y los parámetros de las curvas 1-D). 
Este enfoque teórico-analítico mediante análisis modal no lineal se extenderá, en la Parte 
II de este trabajo, para tener en cuenta la anisotropía desigual y la heterogeneidad del 
miocardio. 


Palabras clave: Modelos de bidominio del miocardio. Diseño de electrodos de marcapa- 
so. Umbrales de excitabilidad. Análisis modal no lineal. Familias de estados umbral. 
Regiones liminales. Estimulación catódica de cierre. Curvas intensidad-duración. 


Rev. SCP” Volumen 18 N*2 2013 pp. 145-174 ISSN 0379-9123 


Abstract: In this first part, after a summary and reformulation of bidomain theory of 
electrocardiology, a simplified mathematical model of the excitation of myocardium by 
pacing electrodes is applied to study the just-threshold dynamics of the electrode- 
myocardium system. A region of influence of the pacemaker electrode on the ventricular 
myocardium is introduced to pose a mixed homogeneous Dirichlet- non homogeneous 
Neumann problem from the nonlinear bidomain field equations, in order to describe the 
polarizing effect of the electrode on the excitable tissue. A nonlinear modal analysis ap- 
proach is used to transform the field equations in a system of nonlinear ordinary differen- 
tlal equations in mode amplitudes. The threshold value of a perturbation that causes insta- 
bility is estimated by a modified version of Eckhaus method. A whole family of thresh- 
old space patterns of membrane polarization is identified in the state space constructed 
from the modal amplitudes. The sizes and shapes of liminal regions (that must be depo- 
larized over threshold of uniform membrane excitation in order to start a propagated ac- 
tion potential) are obtained. Analytical expressions of cardiac strength-duration curves are 
derived, for cathodic controlled current stimulating pulses, including closed form formu- 
lae for time constants (or chronaxies) and rheobases in terms of geometric and electro- 
chemical electrode parameters, electrode-myocardium distance, electric properties of the 
volume conductor of the tissues, and the local distribution of electrotonic and excitability 
parameters of myocardium. Several experimental results are thus explained (liminal re- 
glons and the behaviour of the parameters of S-D curves) and some new facts are predict- 
ed (related with liminal regions, the family of electrode-myocardium threshold states and 
the parameters of the S-D curves). In Part II of this work, the theoretical- analytical ap- 
proach using nonlinear modal analysis will be extended to take into account the unequal 
anisotropy and heterogeneity of myocardium. 


Key words: Bidomain models of the heart. Pacemaker electrodes design. Excitability 
thresholds. Nonlinear modal analysis. Families of threshold states. Liminal regions. 
Cathodic make stimulation. Strength-duration curves. 


(1)Introducción 


Una de las innovaciones médicas más importantes del siglo pasado fue el marca- 
paso cardíaco implantable, introducido hace más de medio siglo para tratar tras- 
tornos del ritmo cardíaco [6] [8] [9] [11] [20]. El desarrollo de los marcapasos 
durante los primeros 40 años se llevó a cabo sin una comprensión a fondo de los 
mecanismos biofísicos subyacentes a la emergencia de un potencial de acción en 
el miocardio. Pese a ello, tanto los ingenieros como los médicos tuvieron bastante 
éxito en definir y concretar las necesarias mejoras en el hardware y en el softwa- 
re, muchas veces por ensayo y error [15], hasta lograr implantar en el paciente un 
verdadero laboratorio electrofisiológico [10]. 

Sin desmerecer los notables resultados alcanzados, hay que reconocer que los 
marcapasos implantables todavía presentan ciertas limitaciones. Por lo general la 
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contracción del miocardio provocada por un estímulo proveniente de un marcapa- 
so externo (a través de electrodos implantados sobre el músculo cardíaco) es me- 
nos eficiente que la contracción generada a partir del marcapaso natural ubicado 
en el nodo sino-auricular, cuando las vías de conducción del impulso funcionan 
normalmente. Dependiendo de las circunstancias del paciente, una contracción 
menos eficiente podría apresurar la instalación de una insuficiencia cardíaca. Este 
es solo un ejemplo de los muchos problemas que se encuentran abiertos, y que en 
su mayor parte involucran de una u otra forma la interfase entre la máquina y el 
corazón: los electrodos de marcapaso [3]. 

Actualmente parece haber un cierto consenso entre los ingenieros y los médicos 
sobre lo siguiente: para avanzar en el diseño, construcción y empleo de los mar- 
capasos externos se requiere una comprensión más profunda de los mecanismos 
subyacentes a la emergencia y propagación de los potenciales de acción en el 
miocardio. En particular, un mejor conocimiento de cómo estos mecanismos bio- 
físicos se relacionan con las arritmias. 

El abordaje reduccionista desde la biología molecular ha conducido a una teoría 
basada en la denominada patología de los canales de membrana, extremadamente 
compleja [32]. Ha permitido plantear tratamientos farmacológicos basados en las 
siempre anheladas “balas de oro”, amparados por una pléyade de patentes de 
invención solicitadas por la industria farmacéutica. Pero los resultados obtenidos 
no parecen acordes con el esfuerzo realizado [32]: es necesario complementar los 
conocimientos sobre la biología y patología molecular de los canales iónicos con 
un abordaje integrador, sistémico, que tenga por objetivo último el intento de 
abarcar el corazón como un todo. 

El abordaje sistémico de la génesis y la propagación del potencial de acción desde 
la fisiología y la fisiopatología ha conducido al desarrollo de varios conceptos 
fundados en experimentos in vivo como el de región liminal de Lindemans [16] y 
a modelos matemáticos cada vez más realistas [6] [32]. Los resultados obtenidos, 
si bien pueden considerarse prometedores, no son todo lo conocidos que cabría 
esperar entre la comunidad de ingenieros y médicos involucrados en el diseño y 
aplicación de marcapasos implantables. 

El propósito principal de esta serie de trabajos es contribuir a una mejor compren- 
sión de los problemas del umbral de estimulación eléctrica del miocardio por 
electrodos externos, a partir del análisis de un modelo matemático del sistema 
electrodo-miocardio. Los resultados de investigación aquí presentados desarro- 
llan, amplían y profundizan las ideas presentadas en la referencia [28]. 

En esta primera parte introducimos y reformulamos los modelos de bidominio de 
modo que resulten aplicables en un análisis de la cardio-estimulación eléctrica 
crónica, sin recurrir a métodos numéricos ab-initio. 
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Luego analizamos los fenómenos de umbral en un modelo simple de sincicio 
excitable no lineal, uniforme e isótropo, con el fin de presentar las herramientas 
analíticas y obtener resultados sobre la excitación del miocardio en un contexto 
libre de las complicaciones que acarrean la no uniformidad y la anisotropía del 
músculo cardíaco. 

Algunas de las características de la excitación catódica de cierre para este modelo 
de sincicio funcional se deducirán y se discutirán, en forma análoga a la empleada 
en las referencias [29] y [30] aplicando herramientas de análisis modal no lineal, 
para analizar los fenómenos de umbral en el caso de una fibra nerviosa excitada 
por un electrodo externo. 

En la segunda parte consideremos las complicaciones adicionales que se derivan 
de la anisotropía desigual y la no uniformidad del miocardio. Estas propiedades 
son fundamentales para poder comprender los mecanismos que subyacen a los 
fenómenos del umbral de excitabilidad eléctrica del corazón, cuando se las consi- 
dera en conjunto con las características específicas de los canales iónicos en las 
membranas excitables del miocardio. En particular, el mecanismo de la excita- 
ción anódica de apertura del miocardio es diferente al mecanismo de la excitación 
anódica de apertura en las fibras nerviosas y en el modelo de sincicio considerado 
en esta primera parte. Debido a las propiedades de los canales de las células car- 
díacas, la excitación anódica de apertura depende de la anisotropía desigual, que 
produce y permite la interacción de cátodos y ánodos virtuales a partir del final 
del pulso estimulador, interacción de la que resulta la excitación [13] [14] [31]. 
En un sincicio isótropo estos electrodos virtuales no pueden producirse, de modo 
que si presenta excitación anódica de apertura, el mecanismo no puede ser el que 
se observa en el caso del miocardio. 

No obstante, como el mecanismo de la excitación catódica de cierre se parece 
más al de las fibras nerviosas, y el efecto de la anisotropía no es decisivo, cabe 
esperar que algunos de los resultados sobre el umbral obtenidos a partir del mode- 
lo simple de sincicio resulten aplicables a la modalidad de estimulación utilizada 
en los marcapasos implantables. 


(2) Modelos de bidominio. Reformulación para el análisis del sistema 
electrodo-miocardio. 


El modelo de cable unidimensional no lineal puede utilizarse para estudiar la 
dinámica del umbral en el caso de fibras nerviosas y musculares estimuladas me- 
diante electrodos externos, debido a la relativa independencia eléctrica que cada 
una de estas fibras presenta respecto de las fibras que la rodean [19]. 

Pero si el tejido estimulado es miocardio o músculo liso del tipo que forma un 
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sincicio funcional, las conexiones eléctricas entre células adyacentes suelen ser lo 
suficientemente fuertes como para que esa relativa independencia se pierda defi- 
nitivamente. En este caso las membranas de las células adyacentes entran en es- 
trecho contacto en varios puntos, compartiendo canales iónicos de gran diámetro 
(conexones) que relacionan entre sí, en forma directa, los espacios intracelulares 
de células miocárdicas vecinas, no solo desde el punto de vista eléctrico, sino 
también desde el punto de vista metabólico [32]. Mientras que el potencial de 
acción que se propaga por una fibra nerviosa puede influir sobre el potencial 
transmembrana de las fibras adyacentes, pero por sí mismo no las excita, la des- 
polarización de una célula de un sincicio eléctrico se propaga a las células vecinas 
de modo tal que si se genera un potencial de acción, éste tiende a transmitirse por 
toda la masa celular eléctricamente acoplada. Para superar esta dificultad se han 
introducido dos tipos de aproximación teórica al estudio de la excitación y la 
propagación del potencial de acción en sincicios funcionales: los modelos de 
autómata [21] y los modelos de bidominio [6] [32]. 

Los modelos de autómata describen los tejidos excitables considerándolos como 
un colectivo de unidades básicas (células) las cuales pueden encontrarse en cada 
instante discreto en uno de un cierto conjunto de estados posibles (por ejemplo: 
en reposo, en estado excitado, en estado refractario). Cada célula puede sufrir un 
cambio de estado al pasar de un instante discreto al siguiente, dependiendo del 
estado de las células vecinas. 

En este trabajo se utilizan modelos de bidominio para describir la génesis del 
potencial de acción en el músculo ventricular (o auricular) en masa [13] [14]. 
Estos modelos describen a tiempo continuo la polarización de los tejidos excita- 
bles y se caracterizan por la introducción de dos medios continuos superpuestos, 
el dominio intracelular y el dominio intersticial, cada uno de los cuales se extien- 
de sobre sobre la totalidad del volumen ocupado por el tejido. Se introducen dos 
campos escalares, los potenciales eléctricos V; (2,7 ) y Volt, ) que se asignan a 


los continuos intracelular e intersticial, respectivamente. Como es usual f repre- 
senta el instante de tiempo y Y el vector de posición del punto del bidominio 
considerado. Se introducen los campos vectoriales de densidad de corriente eléc- 


trica en los dominios intracelular e intersticial J , (1,7 ) y 0 E (2,7 ) respectiva- 

mente. La densidad de corriente eléctrica total en el bidominio J (e, r ) se expre- 

sa como suma de las densidades de corriente eléctrica en cada dominio: 
JT=J+J, (1) 


La conexión eléctrica entre estos los dominios se efectúa a través de un tercer 
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medio continuo formado por las membranas excitables, que también se supone 
ocupan la totalidad del volumen de tejido considerado. A las membranas del ter- 
cer dominio se les asignan dos campos escalares: el potencial local de la mem- 


brana excitable Y, (t,7)=V (c,F)-V,(t,7) y la densidad local de corriente a 


través de las membranas excitables J a, 17) , siendo W un vector cuyas com- 
ponentes son las variables de activación y de recuperación que caracterizan la 
cinética de los canales iónicos, o si se quiere, los efectos cooperativos (transicio- 
nes de fase fuera del equilibrio termodinámico) que se producen cuando las 
membranas y sus zonas adyacentes son perturbadas por un fenómeno de polariza- 
ción localizado. 

Se introduce un campo escalar Xy (e,F ) para representar la fracción de área de 


membranas excitables por unidad de volumen de bidominio y en cada punto se 
postulan las relaciones: 


VeJ,=-%uJy Qa) Vel =+X%u Ju (2b) 


En las ecuaciones (2) el símbolo V e corresponde al operador divergencia. 

El resto de los tejidos, incluyendo los tejidos donde se encuentra implantado el 
electrodo estimulador, se describen como un único continuo, denominado espacio 
extra-miocárdico, al que le asignan un campo de potencial eléctrico D, (1,7 ) y 


un campo de densidad de corriente eléctrica J, o (e, r ). 


Se añaden condiciones en la frontera entre la superficie del electrodo y el conti- 
nuo que representa a los tejidos adyacentes, y entre este último y el bidominio 
que representa al tejido cardíaco ventricular. La frontera entre el bidominio y el 
espacio extra-miocárdico se representa mediante una superficie regular. Se postu- 
lan dos relaciones de continuidad en cada punto de la frontera y en cada instante: 


V,=0, (3) ño J=ñ9b, (4) 
Con 7103) se representa el campo de vectores unitarios (versores) normales a 


cada uno de los puntos de esa superficie frontera y con e el producto escalar. 

Las vías de conducción especializada (el haz de His en el miocardio ventricular y 
las discutidas, y por lo general no admitidas, vías específicas de conducción auri- 
cular) y los núcleos marcapaso (los núcleos sinusal y aurículo-ventricular), histo- 
lógica y funcionalmente bien diferenciados, pueden modelarse por separado me- 
diante un sistema de fibras o haces unidimensionales de fibras que se ramifican a 
partir de un tramo inicial que sale del núcleo correspondiente. 

Considerados como conductores de volumen, tanto el músculo ventricular como 
el músculo auricular son anisótropos y en general, no uniformes. En el modelo de 
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bidominio, la conductividad en cada punto de cada dominio se puede caracterizar 
por un tensor de conductividad. 


Si G, y G, simbolizan los campos tensoriales simétricos que caracterizan la 


conductividad eléctrica en cada punto de los dominios intracelular e intersticial, 
respectivamente, y si f, y f, representan las fracciones de volumen de los domi- 
nios intracelular e intersticial, respectivamente, de modo que f,+f, =1, se pos- 
tulan las siguientes generalizaciones de la ley de Ohm: 

J,=-f, Ge V/, (Sa) ¿"LG ¿0VY, (5 b) 
En estas ecuaciones, como es usual, V simboliza el operador gradiente. 
Puesto que los tensores de conductividad son simétricos, en cada punto de cada 
dominio hay tres direcciones principales ortogonales tales que la representación 
matricial del tensor en esos ejes es diagonal. Así pues, la conductividad tensorial 
local se puede caracterizar mediante tres coeficientes escalares de conductividad, 
uno por cada dirección principal del tensor de conductividad. Como al menos dos 
de estos coeficientes son diferentes entre sí, tanto el dominio intracelular como el 
intersticial son anisótropos tomados por separado. Las direcciones principales del 
tensor de conductividad son aproximadamente las mismas en un mismo punto 
para ambos dominios superpuestos. Pero como los tensores de conductividad de 
esos dominios no son proporcionales en un mismo punto, el grado de anisotropía 
es diferente en el dominio intracelular respecto del dominio intersticial: la aniso- 
tropía es desigual. 
En particular, el miocardio ventricular presenta una anisotropía desigual muy 
marcada que, junto con su no uniformidad, y con las modificaciones que los gra- 
dos de anisotropía y de no uniformidad presentan con el paso del tiempo, dan 
origen a importantes consecuencias electrofisiológicas y electro-patológicas [6] 
[32]. 
La conductividad del espacio extra-miocárdico se describe mediante otro campo 


tensorial simétrico G, y se postula una generalización de la ley de Ohm para 
conectar la densidad de corriente Jo (t,F ) con el gradiente del potencial escalar 
Dole,7): J,=-G,+VO, (6) 

Si bien las ecuaciones de campo y las condiciones en la frontera que se utilizan en 
la teoría de bidominios se postulan sobre la base de condiciones de plausibilidad, 
es posible fundamentarlas mediante procedimientos de homogenización basados 


en la teoría de perturbaciones [13] o en la teoría de campos promediados [24]. 
Para aplicar métodos analíticos en el estudio de la génesis del potencial de acción 
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en el miocardio es conveniente sustituir los dominios intracelular e intersticial por 
dos continuos que permiten una descripción más directa del proceso de excitación 
por medio de electrodos de marcapaso: el continuo formado por las membranas 
excitables (cuyas variables de estado son los campos Vu (t,7 ) y W(t,F)) y un 
nuevo continuo, que se comporta como un conductor de volumen óhmico, gene- 
ralmente no uniforme y anisótropo, caracterizado por un campo escalar (e, ) y 
por un nuevo tensor simétrico de conductividad definido a partir de las fracciones 
de volumen de los dominios originales y de sus tensores de conductividad [24] 


[25]: G=f. -G,+f.-G, (7) 
El campo e, r ) se define de modo tal que la densidad total de corriente 
j (2, r ) en el bidominio verifique: J=-GevY (8) 
Entonces, de (1), (5 a), (5 b), (7) y (8) se desprende: 

GeVY =f Ge VV,+f.-G,eVY, (9) 


Considerados como operadores (díadas) los tensores de conductividad son estric- 
tamente positivos. Entonces de (7) resulta que, considerado como operador, el 


tensor G es estrictamente positivo y por tanto invertible [12] [22]. 
Si G” es el correspondiente operador inverso (que, como G mismo, es también 
un tensor simétrico): vy =07 elf, 6, eVY, +f, 6, SA) (10) 
De (1), Q a) y Q b) se desprende que V e.J =0 de modo que de esta última 
ecuación y de (8) resulta que el campo Y(t,7 ) verifica la ecuación: 
Ve(Gevw(17))= 0 (11) 
En suma, el bidominio viene ahora caracterizado por las variables de estado 
V. (6,7) y W(t,7) del continuo de las membranas excitables que de aquí en 
adelante denominaremos medio rojo, y la variable de estado (e, ) del nuevo 


continuo que denominaremos medio azul. 
Ahora expresamos la densidad de corriente eléctrica a través de las membranas 
excitables como suma de una densidad de corriente de desplazamiento 


O 3 
Ei rá (2,7) (siendo C,, la capacidad por unidad de área de membrana) y 
É 


una densidad de corriente iónica J,,, lr, , 7), se puede deducir la siguiente ecua- 


ción de evolución para el medio rojo [24]: 
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Lu (e, Vale 7) TS (7) =Ve(Covr,)+ve(Dej) (12) 


Por definición los tensores auxiliares simétricos C (que es, además, estrictamente 


positivo) y [DD vienen dados en función de los tensores de conductividad de los 
dominios intracelular e intersticial por las relaciones: 


Ep D=0 "elf. 3,18) (14) 


En (12) el medio rojo se acopla con el azul a través de la densidad de corriente 
eléctrica en el medio azul, producida por el electrodo de marcapaso. El término 
Ve (D e. J ) es el equivalente para el miocardio de la función activante de Rattay 
de la teoría del cable aplicada a la excitación de fibras nerviosas mediante elec- 
trodos externos [29] [30]. Se ha propuesto denominarla “función activante gene- 
ralizada” [31]. 

Cuando la anisotropía es desigual y la densidad de corriente J generada por un 
electrodo externo al miocardio no es nula, la perturbación del campo Vu (t,7 ) 
producida por la función activante generalizada permite explicar la aparición de 
cátodos y ánodos virtuales en el interior del tejido excitable. Los electrodos vir- 
tuales y sus implicancias se estudiarán en la segunda parte del trabajo. 

La descripción de la dinámica del medio rojo se completa con un modelo mate- 
mático no lineal para la membrana excitable en términos de una especificación de 


cómo depende la densidad total de corrientes iónicas /,, de las variables de esta- 


do de la membrana , 3 17) y de cómo evolucionan en el tiempo las variables de 


activación y de recuperación: 


E AA (15) PF 1) (16) 


ion 
Se añade la condición en la frontera (4). La condición en la frontera (3) se puede 
sustituir por una restricción a la densidad de corriente J, ¿ en el dominio intracelu- 
lar en cada punto de la superficie de separación entre el bidominio miocárdico y 
el espacio extra-miocárdico [24]. Para los fines de la presente investigación, 
debido a que la resistencia de las membranas es muy elevada en comparación con 
la resistencia en el interior celular, supondremos que la componente de la densi- 
dad de corriente intracelular normal a la superficie del bidominio miocárdico 
puede despreciarse: ne.J,=0 

Teniendo en cuenta esta última hipótesis, y que a partir de las ecuaciones 
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Y, (t,F)=V (t,F)-V,(t,7), (8) y (9) es posible demostrar que la densidad de 
corriente en el dominio intracelular verifica J de G, . (qa .J )- Cev Vis 0:58 
desprende la condición en la frontera que utilizaremos de ahora en adelante: 

ño (C.VY, )= 01,3, e(G +3) (17) 
SsiñoJ y Se supone conocida, el campo y (e, r ) se determina a partir de las ecua- 


ciones (8) y (11) junto con la condición en la frontera (4). En otro caso es necesa- 
rio resolver simultáneamente las ecuaciones del medios rojo y azul junto con las 
ecuaciones del espacio extra-miocárdico, añadiendo condiciones en la frontera 
entre el electrodo de marcapaso y el conductor de volumen tisular y en el resto de 
la frontera de este último conductor. 


(3) Modelo matemático de sincicio excitable, homogéneo e isótropo. 


Cuando el sincicio eléctrico es isótropo: G,=G,-1] G,=G 1 G=G-1 


dy O (E 
Siendo 1 el tensor identidad yG =f,-G,+f,-G,. Además: G” = E y 


GA. 6, 7” 5 l (f.-G,—f, :G,) 


Entonces C=C-7, C= D=D-[.,D=-. 


LGA EAT.) 
Suponiendo que el sincicio es homogéneo, las ecuaciones de los medios azul y 
rojo se reducen, respectivamente, a: y? ple, Y ) =0 (18) 
19) E Z 
tu [Cs lr) 7)]=0- 7, (19 a) 


Para analizar la dinámica del umbral supondremos que las variables de activación 
se encuentran relajadas al equilibrio con el potencial de membrana, y que las 
variables de recuperación se pueden resumir en una sola variable W : 
J,, =JV,,W) (19 b) 1 = El) (19 c) 

Para describir la dinámica del umbral en las membranas excitables del miocardio 
ventricular y auricular con solo dos variables de estado V,, y W , se puede utilizar 
el modelo que propusieron Aliev y Panfilov [1]. Este modelo es una modificación 
del modelo clásico de Fitzhugh-Nagumo construida especificamente para ser 
aplicada al miocardio contráctil. El modelo de Aliev-Panfilov completo será utili- 
zado en la segunda parte de este trabajo. 
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En esta primera parte supondremos que la variable de recuperación se mantiene 
congelada en su valor de reposo. En ese caso la fórmula para la densidad de co- 
rriente iónica a través de las membranas excitables es la misma en los modelos de 
Aliev-Panfilov y de Fitzhugh-Nagumo. Se expresa solamente en función de la 
diferencia v entre el potencial de membrana V,, y su valor de reposo V, : 


y. 


ion 


=4(9)==—(P=bw +0v) (0) 
M 

Los parámetros by c son positivos, al igual que la resistencia de la unidad de 

área de membrana R,, en las proximidades de V,, =V,. El polinomio cúbico 

posee tres raíces: v=0 que corresponde al potencial de reposo, v=u que co- 

rresponde al umbral de membrana en condiciones de despolarización uniforme, y 

v= yv, que corresponde a un potencial de estado excitado en ausencia de efectos 

debidos a la variable de recuperación. Es la misma fórmula que utilizamos para 

desarrollar un análisis modal en una ecuación de cable no lineal [30]. 

La ecuación (19 a) se reduce entonces a la siguiente: 


0) E 
Tyco MF) += bw +0 = Ap VOY Q1) 
É 
La constante de tiempo de membrana viene dada por: Tu =Ru"Cuq  Q2) 
R 
La constante de espacio del sincicio eléctrico resulta: Ay =. -C (23) 
Xm 
La condición (17) en la frontera entre el bidominio miocárdico y el espacio extra- 
OV,  0V 1 5 
miocárdico se reduce a: —=—= : li e J ) (24 a) 
ón ón f,-G, 
La condición (4) en la frontera permanece igual: noJ=ñeJ o. (Q4b) 


0) 
Como es usual, el símbolo V*representa al operador laplaciano, y E simboliza 
n 


la derivada en dirección de la normal a la frontera entre el bidominio cardíaco y el 
espacio extra-miocárdico. 


(3.1) Análisis modal no lineal, familias de estados umbral y regiones limina- 
les de Lindemans. 

Los fenómenos de polarización de las membranas excitables del miocardio, debi- 
dos a electrodos externos o producidos espontáneamente en focos ectópicos inte- 
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riores al músculo ventricular, son siempre fenómenos localizados. Entonces es 
posible asumir que la polarización inicial v(0, F ) ocupa una región acotada (2 del 
bidominio que representa al tejido excitable. En ese caso es razonable formular la 
siguiente conjetura: la polarización se puede considera nula en la frontera 0(2 de 


la mencionada región durante todo el proceso que, o bien finaliza con el retorno 
del potencial de membrana a su valor de reposo, o bien finaliza en la emergencia 
de un potencial de acción: vlt,F)=0 para FedQ (25) 
Si consideramos, para comenzar, un bidominio que ocupa todo el espacio, con la 
región (2 en su interior, podemos prescindir de las condiciones en la frontera 
(24) para concentrarnos en las soluciones de la ecuación a derivadas parciales no 
lineales (21) en el interior de (2 y con la condición de borde (25). 


Busquemos construir soluciones de la forma: vt, r ) = y A, (e) 4, (7 ) (26) 


Si la región Q es lo bastante regular, podemos emplear las funciones propias del 
operador —V? con condiciones de Dirichlet homogéneas sobre 062, es decir 


- Vu, = TE -U,, para representar la dependencia espacial de la solución para el 
potencial de membrana [17] [18] [22]. Como la región es tridimensional, en ge- 
neral aparecen tres índices enteros positivos /,m,n que, para simplificar la nota- 
ción, resumimos en un único índice 7». 

Las funciones propias forman una base ortogonal para el espacio de las funciones 
de cuadrado integrable en el sentido de Lebesgue (entre las cuales se encuentra la 
solución vlt,F )) y los valores propios HE son positivos y forman una sucesión 


creciente con límite infinito. De ahora en adelante supondremos que las funciones 
de la base están normalizadas [17] [18]: Ju,,(*)-u,(+)-dV =5,,, 
Q 


En forma análoga a la empleada en la referencia [29], nos a determinar 


una solución aproximada (ansatz) Y (e, 7) 24d) ) de la ecuación (21) 


aplicando un método de Galerkin [17]. 
Se obtienen este sistema de ecuaciones diferenciales ordinarias no lineales para 
las amplitudes de modo 4 (e ) (n=12,...,M ): 


da, 
Ty 2 =-(1+ Ay 4, )A, YO. A,A, -0 y Bavor Ay Ay A, (27) 


pq=1 pqr=1 


Los coeficientes O vienen dados por las siguientes integrales de volumen: 


Rev. SCP” Volumen 18 N*2 2013 pp. 145-174 ISSN 0379-9123 


dE j u,(+)-u,(+)-u,(F)-dV (28) Osa J u,(*)-u,(*)-u,(*)-u,(*).dV (29) 


Al igual que en el caso de la fibra excitable, estos coeficientes dependen solamen- 
te de las funciones propias del operador. En general, siempre pueden calcularse 
numéricamente. Si la región ( es lo bastante simétrica (paralelepípedo rectangu- 
lar, esfera o elipsoide) podrían obtenerse mediante el uso de paquetes de software 
para cálculo simbólico. 

El estado de reposo de las membranas del medio rojo corresponde al origen del 
espacio de las amplitudes de modo. 

Supongamos que las membranas excitables se encuentran en reposo. Un campo 
de corrientes eléctricas en el conductor de volumen de los tejidos, producido por 
un electrodo de marcapaso (y acoplado con el medio rojo a través de las condi- 
ciones en la frontera (24)), o las corrientes locales asociadas a la aparición espon- 
tánea de un foco ectópico producen una perturbación no uniforme que aleja a las 
membranas de su estado de reposo. 

Cuando la amplitud de la perturbación es lo bastante pequeña, cabe esperar que el 
estado del medio rojo retorne directamente al reposo. Esto se comprueba fácil- 
mente aplicando el primer método de estabilidad debido a Poincaré y Liapunov 
junto con el teorema de Hartman-Grobman [4] [17] [18]. (Comenzando por la 
aproximación lineal a las ecuaciones (27) en torno al origen del espacio de las 
amplitudes de modo). Pero, y esto es lo más interesante desde el punto de vista 
de la estimulación eléctrica, cabe esperar que si la amplitud de la perturbación es 
lo bastante grande, el estado del medio rojo se aleje aún más del reposo y se pro- 
duzca la excitación. 

Si la dinámica del medio rojo se describe mediante el sistema (27), la amplitud de 
una perturbación que desestabiliza el estado de reposo puede estudiarse aplicando 
una modificación apropiada de un método desarrollado por Eckhaus [4] [5] [29]. 
Un análisis numérico análogo al realizado en la referencia [29] sugiere que el 
modo de amplitud 4, (s) correspondiente al valor propio mínimo 4, es dominan- 


te. Desacoplando el modo dominante de los modos de orden superior se obtiene: 


dA 
Tu rn dd 41) A, + b-0,, Aj CG Ay (0) 


2081111 Gun 


Si ese discriminante es positivo, la ecuación diferencial (30) posee tres puntos de 
equilibrio: 


2 
Definamos ahora el discriminante A = E) qe li +A ) (31) 
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NA Bb: 


Lu 


JA (32b) 427 + Ya (320) 


2001111 Có1 


A =0 (322) A 


En este caso Á , , es un estado de equilibrio instable. Si la amplitud inicial se en- 


cuentra entre A,=0 (el estado de reposo de las membranas excitables del sinci- 


cio) y A, la amplitud del modo dominante retorna siempre al estado de reposo. 


lu > 
A este intervalo, igual que en la referencia [29], lo denominaremos conjunto de 
decaimiento. 


Si la amplitud inicial se encuentra comprendida entre A, y A ¿la amplitud del 


Lu 
modo dominante tiende siempre hacia el estado excitado A, ,. A este intervalo lo 
denominaremos conjunto amplificador. 

Por tanto, A, es una amplitud umbral que debe ser superada para que la exci- 
tación del sincicio pueda producirse. 

La solución, aproximada hasta el modo dominante es %(t,7)= 4 (t)-u (7). 
Entonces, siempre dentro de la aproximación del modo dominante, el perfil esta- 
cionario 7, (F)= A, 4 (+) puede considerarse como un perfil umbral. 

La Figura 1 muestra un diagrama de bifurcación de las tres amplitudes de equili- 
brio en función de un parámetro proporcional a 4, . 


Figura 1. Diagrama cualitativo que muestra una bifurcación global de las amplitudes de 
equilibrio del modo dominante en función del valor propio del modo dominante. 


La línea punteada vertical en la Figura 1 asciende desde el valor crítico del pará- 
metro P que corresponde a 4 . . 


si v(0,*)>7, (7) la excitación se produce. Como se impuso la restricción (25), 
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el potencial de membrana vlt,F ) se aproxima asintóticamente, para t>00, al 
perfil estacionario %,(7)= 4, -u(7). Si v(0,7)<%,(F) , vlt,F) tiende a cero en 
todos los puntos de la región (2 : la excitación fracasa. 

A medida que Á se aproxima a cero desde valores positivos, las amplitudes um- 


bral y de estado excitado se aproximan, hasta que, para A = 0 coinciden. Cuando 


A se hace negativo desaparece el umbral y el sincicio ya no puede ser excitado. 
Supongamos que los parámetros de la membrana y la constante de espacio del 
sincicio permanecen invariantes y que para valores lo suficientemente próximos a 


2 . oo. “> . .,s 
cero de 44, el discriminante es positivo. Entonces existe un valor critico 1, . 
tal que Á se anula si 4, = 41, ., y se hace negativo cuando 4, >p4,.- 


De (31) se deprende que: 


1 , OO 
1 E 


1 (33) 
Los valores propios del problema — V?u, = rs -4, en (Q, con condiciones ho- 
mogéneas de Dirichlet sobre 042, dependen solo de la región Q.. 


, ; y n:«T , 
Si Q es una esfera de radio 7, los valores propios son 4, =—— y las funciones 
/ 
0 


propias verifican [17] [18]: ro) 22] 


h 
Si Qes un paralelepípedo rectangular de lados 2-a, 2-by 2-c los valores 


2 Ze 2 2, 
y 2» T l mon . ; y 
propios SOd 4. ==" AR mientras que las funciones propias 


4£la be 


L-T-x m-T-x N-T-X 
17] [18]: =Si -si -si 
son [17] [18] Uy mn (2%, y, 2) sin E ) sin 2.5 ) sin 20 ) 


Si Qes un elipsoide, el problema se puede resolver en términos de funciones de 
Lamé separando variables mediante coordenadas elípticas. Si (. es un cilindro 
aparecen productos de funciones de Bessel y trigonométricas, cuando se utilizan 
coordenadas cilíndricas para separar variables. Pero en general, para una región 
de forma arbitraria no es posible reducir el problema a uno que involucre algunas 
de las funciones clásicas de la fisica-matemática. No obstante, es posible acotar 
inferior y superiormente los valores propios recurriendo a regiones para las cuales 
el problema posee una solución analítica conocida, sobre la base que suministra 
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un teorema que aparece expuesto en las referencias [22] y [25]. En suma, ocurre 
lo siguiente. Supongamos una familia de regiones Q(s), donde € es un paráme- 
tro real (o un conjunto de tales parámetros). Supongamos que resolvemos el pro- 
blema -V?u(r) =p -u(r) sobre O(s) con la restricción u(F)= O para todo 
r€ 0Qls), obteniendo un conjunto de funciones propias u, (7 2) y de valores 
propios correspondientes 1, (5) , 

Si las dimensiones de las regiones dependen en forma regular del parámetro e, lo 
mismo ocurre con los valores propios 1, (e) [22]. 

Entonces es posible demostrar que si O(s,)= Ole) Qls,), donde la inclusión 
es en sentido estricto, los valores propios correspondientes al mismo índice n (o 
índice múltiple /,m,n ) verifican la desigualdad estricta yu *(s,)> un, (s)> u,(e,). 


Por ejemplo, si Q(s) es un elipsoide de semiejes a, by c podemos tomar como 


2-4 
región Q(s,) el paralelepípedo rectangular con el mismo centro y lados ——- 


a 


2b 26 ¿ 
Ed 5 y como región Q(s,) el paralelepíipedo rectangular análogo pero 


con lados 2-a, 2-by 2-c. Entonces el valor propio HimnlE) del problema 


planteado sobre el elipsoide verifica Mom (Es) S Mim (E) S Mim.» (€), 


, HEN AR 2 SR AS 
es decir: a : + + <Himn (e)< 1 ES 
cm a 


2 2 2 2 
a Doe 4 be € 
Como los índices son enteros positivos, el valor propio del modo dominante es 


2 
A le)s 1 (s) y queda acotado inferiormente por y ( s + z + j ) y supe- 


O O 
y á 31 (1 s 1 :? 1 
riormente por . , 
4 ae De 
Con suficiente aproximación para los propósitos de este trabajo: 
2 
2 S:7% 1 1 1 
pan 8 (Fr) 5 


Para que una polarización inicial de las membranas del elipsoide pueda originar 


$ so SER 2 2 : 
una excitación del sincicio, /e;, , debe ser menor que 4”, . Como consecuencia 


las dimensiones del elipsoide deben verificar la relación: 
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2 
1 1 1 8 ; Me 
OS 2 

a yO Sa 
Consideremos ahora un caso más interesante desde el punto de vista de la cardio- 
estimulación por medio de marcapasos, en cual la región R ocupada por el sinci- 
cio es acotada y la polarización alcanza a su frontera OR. Una vez finalizado el 


(35) 


pulso estimulador, la frontera O(2 de la región ( ocupada por la polarización 
resultante se puede dividir en dos partes. Una de ellas, LDANOR se encuentra en 
la superficie que separa el bidominio cardíaco del espacio extra-miocárdico. La 
otra, OQ 00 Nn OR)se encuentra en el interior del bidominio cardíaco. 


La solución aproximada v(t,7 )- y A, (2)-u, (7) se determina aplicando un mé- 


n=1 


todo de Galerkin, igual que antes. Pero las funciones propias u, (7 ) que aparecen 


. 2 2 . 
en este ansatz son las soluciones del problema — Wu = u“-u en (2 con condi- 


ciones homogéneas mixtas de Dirichlet y de Neumann en la frontera DA: 
u(r)=0 sobre 0NQÁ(DANAOR) (Dirichlet) y als) =0 sobre 
n 


DOMO R (Neumam). 

La condición homogénea de Neumann se obtiene a partir de la condición en la 

frontera (24 a) para una densidad de corriente eléctrica nula en el medio azul. 

Las ecuaciones (27) de evolución de las amplitudes de modo se aplican también 

en este caso, pero con los valores propios correspondientes al problema mixto. 

Si R fuera, en coordenadas cartesianas otrogonales, el hemi-espacio z<0, su 

frontera OR sería el plano z=0. Supongamos que la región (2 es un hemi- 
2 2 2 


<1l para z=0. En este caso, por simetría, los valores pro- 


2 


: : A o, XK 
pios son los mismos valores propios para la región E + 


2 


2 
z z , 
Le += Sl interior al 
C 
elipsoide completo, contenido en el sincicio no acotado que estudiamos antes. Las 
funciones propias también son las mismas, con su dominio restringido al hemi- 


elipsoide. También se mantiene la desigualdad (34) como condición para que una 
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polarización inicial de las membranas del hemi-elipsoide pueda originar una excl- 
tación en el hemi-espacio ocupado por el sincicio. 


El valor crítico Ha para el hemi-elipsoide es el mismo valor crítico para el elip- 


soide correspondiente, pero el volumen de uno es la mitad del volumen del otro. 
Convengamos en denominar región crítica (2, a una región tal que una despola- 


rización en una región estrictamente incluida en (2, no puede producir la excita- 


ción del tejido, mientras que una despolarización en una región en cual (2. está 

estrictamente incluida, puede producir excitación. El valor propio 1,” corres- 
. ., rn. rus 2 z 

pondiente a una región crítica es el valor crítico 4, dado por la fórmula (33). 


En esa fórmula los coeficientes O dependen de la región (2, a través de las fun- 


ciones propias que aparecen en las integrales de volumen (28) y (29). 
Para cada región de influencia despolarizada por acción de un foco ectópico o de 
un electrodo externo, que incluya estrictamente a una (2, , hay un estado umbral 


del sincicio excitable 5, (7)= 4,,-u (7). 


Hay muchas regiones de influencia en esas condiciones, y a cada una le corres- 
ponde un patrón espacial umbral de despolarización de las membranas excitables. 
Por tanto no hay un único estado umbral, sino una familia de estados umbral. 

Por último, el concepto de región crítica puede aplicarse para interpretar un ha- 
llazgo experimental de Lindemans: la región liminal. Para generar un potencial 
de acción capaz de propagarse es necesario despolarizar por encima del umbral 
uniforme de membrana un volumen o masa crítica de células cardíacas [16]. 


(3.2) Estimación de la forma y dimensiones de la región de influencia de un 
electrodo de marcapaso sobre el miocardio. 

Consideremos, para concretar un electrodo bipolar situado a cierta distancia del 
miocardio. El músculo cardíaco se representa por un sincicio homogéneo e 1só- 
tropo que ocupa en hemi-espacio z>0. 


La Figura 2 muestra el electrodo bipolar con el cátodo proximal al tejido miocár- 
dico. El miocardio excitable se representa por un hemi-espacio ocupado por el 
medio rojo (las membranas excitables) superpuesto con el medio azul (el conduc- 
tor de volumen equivalente al tejido cardíaco). 

Se observa un anillo anódico distal respecto del bidominio, y un hemi-esferoide 
que simboliza la región de influencia del electrodo bipolar sobre el miocardio. Se 
representan algunas líneas de corriente. Todas ellas que salen del anillo anódico 
para terminar en el cátodo. 
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Anillo anódico xx 


; E ———4 5 << > 
Líneas de —— ; 
corriente 234 ES 
eléctrica _ pr” —— 


< Región de >) 
HS R «7 influencia / 


Sincicio eléctrico 


Figura 2. Región de influencia de un electrodo de marcapaso sobre un sincicio eléctrico 
homogéneo e isótropo: un hemi-esferoide de radio R, y semieje vertical c, . 


El campo de corrientes eléctricas producido por el electrodo se acopla con la 
ecuación (21) que describe la polarización de las membranas excitables, a través 
de las condiciones en la frontera (24 a) y (24 b). En la frontera z=0 del hemi- 
espacio la componente normal de la densidad de corriente eléctrica en el conduc- 
tor de volumen se puede representar así: 
ño J=ñeJ, =F(po)-i(t) (36) 

Debido a la simetría asumida para el campo de corrientes producido por el elec- 
trodo bipolar, el factor geométrico F (p) depende solamente de la distancia 
p del punto considerado del plano z=0al origen de este plano. Debido a que 
los campos se pueden considerar cuasi-estacionarios, la corriente eléctrica global 
introducida a través de electrodo aparece como un factor en (36). Entonces, de 
las ecuaciones (24) se obtiene: 

Ov 1 5 

— ==—— (ne J, Ju Flo): il: (37 

gero) ) 
La Figura 3 muestra un esquema cualitativo del factor geométrico. La función 
y=F(p) es simétrica y disminuye aproximándose a 0, sin cambiar de signo, a 
medida que la distancia del punto considerado al origen tiende a infinito. 
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Figura 3. Esquema de la dependencia radial de la componente de la densidad de corriente 
paralela al eje z en la superficie z=0 del sincicio. 


Es posible proceder como se hizo en las referencias [29] y [30] para el caso de 
una fibra excitable, introduciendo ahora un alcance radial de la despolarización 
debida al electrodo. A partir de la intersección de la recta tangente a la curva 


y=F (p) en el punto de inflexión / de abscisa p,, con el eje de abscisas (ver la 
Figura 3) se puede definir un radio R, a partir del cual convenimos en despreciar 
el efecto del electrodo sobre las membranas del sincicio eléctrico: 


E(o ) 


(38) 


El radio (38) es el radio del esferoide de influencia del electrodo. Por ejemplo, si 
el electrodo fuera un cátodo puntual, situado sobre el eje Z a una distancia Z, del 


Zo 


sincicio, y con el ánodo en el infinito: F(p)= Para este 


3 
47 (+2 


caso p =p, = 5 El radio del esferoide R, = p;¡+ =—"Z¿y  €S pro- 


porcional a la distancia del electrodo al tejido excitable. 

Asumiendo que la región de influencia es un esferoide, solo resta estimar el semi- 
eje C, para definir sus dimensiones. 

Teniendo en cuenta la simetría, la ecuación (21) se puede re-escribir así, em- 
pleando coordenadas cilíndricas p,Z: 
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O v Il 2 "ia NS 
A E (p HKÑ (39) 
O t pop Op. 0Z 
Supongamos que v(0, p,z)=0 (es decir, el estado inicial de las membranas del 
medio rojo es el estado de reposo). Entonces, teniendo en cuenta el radio del esfe- 
roide de influencia, cuando el electrodo comienza a actuar a partir del instante 
inicial, podemos suponer que el potencial de membrana se anula para p=R),: 
v(t,R,,z)=0para todo t>0 y todo z>0. 


Buscamos una solución de la forma: 


vít, p,2)= 9 An(t,z) v,(0) (40) 
n=1 
Las funciones propias v, (/p) son solución del problema de Sturm-Louwville: 
ld dv 
(0 (p)=0, v, (0) (41) Y. (Rg)=0 (42) 
pdp” dp 
Las soluciones se pueden expresar así: v,(p)=C, J,.(O, p) (43) 


J,(G) es la function de Bessel de primera especie y orden cero [17] [18]. 


Como v,(R,) = 0, a partir de (43) se deduce: o. => (44) 


Aquí 4, <A , <A ¿ <... simbolizan los ceros de esta función de Bessel. 
Los cuadrados de los primeros tres ceros de las funciones de Bessel son: 
a, =5.7830  0a4,=30.4715 0, =74.8865 
Debido a la hipótesis según la cual el primer modo es el dominante, cabe esperar 
entonces que la máxima profundidad de penetración se asocie a la función propia 
de orden más bajo v,(pP). Además cabe suponer que se alcanza cuando la polari- 
zación de las membranas corresponde a un estado estacionario sub-umbral, debi- 
do a un escalón de corriente polarizante aplicado por el electrodo bipolar. 
Sustituyendo la solución tentativa (ansatz) estacionaria A,.,(z)-v,(p) en la 
ecuación (39) y tomando la aproximación lineal a la expresión resultante, obte- 
nemos la siguiente ecuación aproximada para A,, (z): 
2 A 
Mt A A (o) (45) 


M 
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La solución con significado físico de esta ecuación debe tender a cero cuando 


VA 


z>-00, por lo cual A,,,(z)=A,,,(0) en, 


y) 
Por definición: L,= = > (46) 
yI+2y 0, 
z 
5.7830 
De (44) se desprende que: 142) % =2 5 = A (47) 
R, R, 


El parámetro L, se puede tomar como una medida de la profundidad máxima de 


penetración de la polarización de las membranas excitables del sincicio (es decir, 
en el medio rojo), es decir, como una medida de la profundidad de influencia del 
electrodo. Teniendo este último en cuenta, introducimos la siguiente longitud 


para el semieje: =P L, (48) 
En (48), fP' es un factor de proporcionalidad que en primera aproximación puede 
tomarse igual a la unidad. 
La Figura 4 muestra en forma cualitativa A como función de ye ; 

M M 


A medida que el radio de influencia crece, L, aumenta aproximándose a A y. 


0 R 


Figura 4. Esquema de la relación entre la profundidad de influencia de un electrodo en 
función del radio de influencia sobre la frontera plana del sincicio eléctrico. 
La fórmula (34) para el valor propio del modo dominante se reduce, para el caso 


z 4 : 2 A S-T a 2 1 
de un hemi-esferoide de influencia a: 4, = al Does dar 
8 Ro C 
Teniendo en cuenta (46), (47) y (48): 


E, 2 ,1f_1 57830 , 
1 8. REAL RA (49) 
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Entonces la condición (35), que debe verificar una región de influencia para in- 
cluir estrictamente a una región crítica, y por tanto para que el tejido pueda ser 


2. 1( 1 57830) 8:14? 
2 de 2 2 Ss 2 < 2 
Ro POLA R, S-T 


Esta desigualdad define un radio de influencia crítico R, . por debajo del cual la 


excitado, se reduce a: 


(50) 


excitación de las membranas no va a ser posible. 


(4) Deducción de fórmulas analíticas correspondientes a las curvas 
intensidad-duración para la estimulación catódica de cierre. 


En (3.1) hallamos las ecuaciones de evolución en el espacio de las amplitudes de 
modo para un sincicio que se supuso era no acotado. En esta parte extenderemos 
ese análisis al caso de un sincicio homogéneo e isótropo acotado por una frontera 
que lo separa del espacio extra-miocárdico, espacio donde se encuentra el elec- 
trodo estimulador (por ejemplo, como el que muestra la Figura 2). Igual que en 
(3.1), representamos la región de influencia del electrodo mediante el símbolo por 
Q y su frontera mediante OO. Esta frontera se puede dividir en dos superficies 


regulares: 02, que se encuentra en el interior del tejido miocárdico (la superfi- 
cie curva del hemi-esferoide) y 4 (2, que se encuentra en la frontera entre el mio- 
cardio y el espacio extra-miocárdico (la superficie circular del hemi-esferoide). 
En los puntos de 00, asumimos que v(t,7)=0. 


Como suponemos que los campos de corriente son cuasi-estacionarios, la densi- 
dad de corriente en el medio azul del bidominio se puede escribir así, introdu- 


ciendo un campo vectorial F (1) que describe la distribución espacial de la co- 
rriente eléctrica global 1(t): J (t, 1) =1(0)- F(5) (51) 
De la condición en la frontera del sincicio (24a) y de la ecuación (51) resulta: 


19) ve.) _ ig) ¡A FE) 


(52 

On f. G. 

Una vez que el pulso estimulador ha finalizado, de (52) se desprende que 

O vt,r j ; a 

Zane =0 sobre 0 (2,. Para obtener las ecuaciones de evolución en términos 
n 


de las amplitudes de modo, utilizamos nuevamente el  ansatz 


(aa pero ahora supondremos que las funciones u, (+) del 
n=1 
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operador verifican una condición homogénea de Dirichlet sobre 4'(2, y una 
condición homogéneas de Neumann sobre/ (2,. Teniendo en cuenta que las 


funciones propias se encuentran HE y son EL y la identidad de 


Green Js V gdV= le V"fdv+ ne ft —dS- 16 , Sustituimos el 


ansatz en la ecuación 20 y luego E varias e obtenemos el si- 
guiente sistema no lineal y no homogéneo de ecuaciones diferenciales ordinarias: 


dA 
Ty EIA A Hb O AA e OA SAGA, 
dl dt d ds E zo ds (53) 
+Ay of, cit) 
; O vit,r 
El término fuente viene dado por: Ft) = J y, PAT) y (54) 


ón 


La magnitud del coeficiente f, mide la intensidad de la perturbación que produce 


00, 


el electrodo estimulador sobre la enésima amplitud de modo. De (52) y (54) se 
obtiene una fórmula que den calcular f,: 


== Jul (+). (ñ, e F(7)).ds (55) 
La G, 00, 

Si la superficie del cátodo es convexa, o si se halla lo bastante alejado del mio- 
cardio como para que las líneas de corriente se abran al aproximarse a la superfi- 
cie del sincicio, cabe esperar que la amplitud del primer modo sea dominante 
[25]. Desacoplándola como se hizo en (3.1) se obtiene: 
Tu (14 A AY A, +b:0,, Aj = 0:05, Aj +4 A (1) (56) 
Cuando finaliza el pulso catódico se anula el término forzante y la dinámica de la 
amplitud del modo dominante satisface la ecuación (30), cuyas propiedades fue- 
ron estudiadas en (3.1). Linealizando (56) en el conjunto de decaimiento, donde 


2 

Aj(t)< A, ,, se obtiene: A «A, y 2 + f, 10) (57) 
dt t. E 

En esta ecuación aparece una constante de tiempo:  = A (58) 


14 E ae 
Integrando entre £=0 y £=d (la duración del pulso catódico), y suponiendo que 
las membranas excitables parten del estado de reposo A,(0)=0, se deduce la 
siguiente fórmula para la amplitud del modo dominante: 
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M Ss 


d A 
A (d)= J Gld —1)-i(t).dt (59) con a) Bl. pe E (60) 


De acuerdo con lo visto en (3.1) un pulso de duración d será umbral siempre y 
cuando A, (d )= Ay, - Aplicando esta última condición a un pulso rectangular 


umbral de amplitud Í, (d ) se obtiene la fórmula de Lapicque-Hill-Blair para el 


I 
miocardio: I, (d ) = al (61) 
1—exp| -— 
+A yn 1) A 
La reobase de corriente verifica: AS M se 4. (62) 
A M EN 
De (49) y (58), con P=1 se deduce la siguiente expresión de la constante de 
Ty 
tiempo: E 2 
a+b- Es 
r, +e 


(63) 
Para deducir la fórmula (63) se utilizó la fórmula Ry =Y* (7 + e) para estimar el 


radio de influencia Ry debido a un cátodo convexo proximal al tejido miocárdi- 
co. Ese radio resulta proporcional a la suma de un radio representativo del elec- 
trodo ?, (el denominado radio equivalente [25] [27]) con el espesor e del cas- 


quete de tejidos no excitables interpuesto entre el electrodo y el miocardio. 
Durante el proceso inflamatorio que se produce en el miocardio adyacente al 
electrodo, desencadenado por la implantación del electrodo, se forma un casquete 
de tejidos no excitables de conductividad eléctrica y espesor variables. 

El espesor de € este casquete reactivo primero crece y luego decrece hasta alcan- 
zar un espesor estable. Comparado con la conductividad eléctrica del miocardio 
próximo al electrodo, la conductividad del casquete también aumenta en la prime- 
ra etapa de la reacción inflamatoria, y luego disminuye cuando se forma la red de 
fibrina. El factor de proporcionalidad Y varía próximo a 4/3 para un cátodo con el 


ánodo en el infinito, aumentando durante la primera fase de la reacción inflamato- 
ria y disminuyendo luego. Ese mismo factor varía en forma análoga próximo a un 
valor algo mayor en el caso de un electrodo bipolar. Los parámetros sin dimen- 


siones 4 y D toman valores próximos a 6 y 38, respectivamente. 
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De (58) y (62), para [4 =1 y teniendo en cuenta que f, oc FR , se obtiene para 
0 


la reobase: LR ( Ed ey E3y0 ql Asu lr =E e) (64) 


u, 
(5) Discusión y Conclusiones 


-El congelamiento de la variable de activación en el caso de las membranas exci- 
tables del miocardio contráctil ventricular y auricular, para el estudio de los um- 
brales de excitación catódica de cierre es una hipótesis menos restrictiva que 
cuando se aplica a las membranas del nervio. Luego de una pequeña recuperación 
al final del frente del potencial de acción, aparece la meseta de despolarización 
que se mantiene durante al menos varias decenas de milisegundos, o sea un orden 
de magnitud por encima de la duración de los pulsos catódicos de los marcapasos 
implantables comunes [2] [19]. 

-La excitación anódica de apertura y la excitación catódica de apertura, que se 
observan en el caso del miocardio [6], no pueden ser estudiadas con el modelo de 
sincicio homogéneo e isótropo construido en esta primera parte. No basta con 
añadir una variable de recuperación al modelo de membrana, como se hizo para 
el análisis modal del umbral de excitación de fibras nerviosas en la referencia 
[30]. 

-La fórmula (63) para la constante de tiempo de la curva intensidad-duración por 
excitación catódica de cierre explica la mayor parte de los resultados experimen- 
tales disponibles, y predice algunos resultados nuevos que tendrían que ser con- 
trastados contra los datos obtenidos en experimentos específicos [25] [26]. Puede 
utilizarse durante el proceso de diseño de electrodos destinados a la cardio- 
estimulación eléctrica crónica. 


- Si Z, =T,+€ no es demasiado pequeño, entonces, como sugiere la Figura 1, la 
amplitud umbral varía poco con Z,. En ese caso la relación (64) predice que la 
reobase aumenta aproximadamente como el cuadrado de z,, pero de tal forma 
que el valor extrapolado de la reobase para z, Y 0 no es nulo. Estos resultados 


coinciden con lo que se obtiene a partir de los datos experimentales disponibles 
[3] 20] 24] [25]. 

-En la segunda parte de este trabajo consideraremos cómo se modifica las expre- 
siones para t, y para 1, ¿ cuando se tiene en cuenta la anisotropía del miocardio. 


-La fórmula (59) para la evolución temporal de la amplitud del modo dominante 
junto con la condición de umbral A, (d )= A, , no son otra cosa que una concre- 


Lu 


ción, para el miocardio considerado como un bidominio, de la teoría clásica de la 
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excitación debida a Blair [25]. Pero a diferencia de lo que acontece con la teoría 
fenomenológica de Blair, la reobase y la constante de tiempo se obtuvieron aquí, 
a través de las fórmulas (58) y (62), como funciones bien definidas de parámetros 
electrofisiológicos y geométricos del sistema electrodo-miocardio. En la segunda 
parte se verá cómo todo esto puede extenderse para tener en cuenta la anisotropía 
y no uniformidad del miocardio. 

-El resultado obtenido para las regiones críticas elipsoidales sugiere la siguiente 
conjetura: el tejido excitable es más vulnerable frente a un foco ectópico que des- 
polariza las membranas de un sincicio acotado en una región adyacente a la fron- 
tera, que frente a un foco ectópico de la misma intensidad que despolariza las 
membranas en una región alejada de la frontera de igual volumen. En la segunda 
parte estudiaremos el efecto de la anisotropía desigual y la no uniformidad del 
miocardio sobre la forma y dimensiones de las regiones críticas, retomando luego 
la vulnerabilidad del miocardio respecto de posibles focos ectópicos. 

-El análisis modal no lineal se pudo llevar a cabo debido a la introducción del 
concepto de región de influencia de un electrodo externo sobre el miocardio. Una 
consecuencia del análisis modal es el concepto de familia de estados umbral. 
Para un mismo corazón, en un mismo estado fisiológico o fisiopatológico, 
la familia de estados umbral que se aplique depende de la forma, dimen- 
siones y distancia del electrodo al miocardio. Desde este punto de vista 
puede distinguirse un concepto global del umbral de excitabilidad del mio- 
cardio de un concepto local de ese umbral. El umbral global depende del 
sistema electrodo-miocardio considerado como un todo (incluyendo el 
estado del casquete de tejidos no excitables interpuesto) y se determina a 
través de las curvas intensidad-duración obtenidas al implantar el electrodo 
o ya en una fase crónica de la estimulación [7] [8] [23]. 

El umbral local es el que se puede determinar a través de esas mismas cur- 
vas, pero obtenidas (por ejemplo) a partir de experimentos de patch-clamp 
con trozos de membrana uniformemente polarizados [32]. 
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eléctrica del miocardio mediante electrodos de marcapaso. Parte II: 
Sincicio eléctrico anisótropo y heterogéneo. 
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Resumen: Esta segunda parte comienza con la reformulación de la teoría de bidominio de la elec- 
trocardiología planteada en el trabajo previo. Se extiende el enfoque teórico-analítico mediante 
análisis modal no lineal para tener en cuenta la anisotropía desigual y la heterogeneidad del miocar- 
dio. Se tienen en cuenta las variables de recuperación de las membranas excitables del miocardio 
empleando un modelo matemático simplificado de la excitación. Se estudia nuevamente, justo hasta 
el umbral y teniendo en cuenta la anisotropía desigual, la dinámica del sistema electrodo-miocardio 
estimulado a través de electrodos de marcapaso. Se re-examina el concepto de región de influencia 
del electrodo de marcapaso sobre el miocardio anisótropo y se tiene en cuenta la función activante 
generalizada para construir la región de influencia. Luego se plantea un problema mixto, homogé- 
neo de Dirichlet y no homogéneo de Neumann a partir de las ecuaciones no lineales del modelo de 
bidominio, de modo de describir el efecto polarizante del electrodo, incluyendo la aparición de 
ánodos y cátodos virtuales. Al igual que en la primera parte, se emplea un enfoque basado en el 
análisis modal no lineal para transformar las ecuaciones de campo en un sistema de ecuaciones 
diferenciales ordinarias no lineales en las amplitudes de modo y el valor umbral de la perturbación 
que causa inestabilidad se estima empleando una versión modificada del método de Eckhaus. Para 
calcular las funciones y los valores propios del operador de Sturm-Louville que permite definir las 
amplitudes de modo se utiliza un espacio virtual (alibi) en el cual el miocardio ventricular puede 
representarse en forma aproximada mediante un modelo de sincicio eléctrico homogéneo e isótropo. 
Se identifican patrones espaciales de polarización de las membranas que corresponden a una fami- 
lia de estados umbral en el espacio de estados construido a partir de las amplitudes de modo, en el 
marco del modelo generalizado y se obtienen las formas y los tamaños de las regiones liminales 
(que deben ser despolarizadas por encima del umbral uniforme de membrana para que pueda gene- 
rarse un potencial de acción propagado) teniendo en cuenta la anisotropía del músculo cardíaco. Se 
analizan las estimulaciones catódica y anódica de cierre y de apertura. Se deducen fórmulas para las 
curvas intensidad-duración correspondientes a la estimulación catódica de cierre, incluyendo fórmu- 
las analíticas para las constantes de tiempo (cronaxias) y reobases, en términos de parámetros geo- 
métricos y electroquímicos del electrodo, de la distancia entre el electrodo y el miocardio ventricu- 
lar, de las propiedades eléctricas del conductor de volumen de los tejidos, y de la distribución local 
de los parámetros electrotónicos y de excitabilidad del músculo cardíaco. Se explican varios resul- 
tados experimentales, se confirman algunos de los resultados de las simulaciones digitales de la 
dinámica del umbral basados en la discretización de las ecuaciones de bidominio y se predicen 
algunos fenómenos nuevos. 


Palabras clave: Modelos de bidominio del miocardio. Diseño de electrodos de marcapaso. Umbra- 
les de excitabilidad. Anisotropía desigual. Heterogeneidad. Alibi. Análisis modal no lineal. Familias 
de estados umbral. Regiones liminales. Ánodos y cátodos virtuales. Estimulación catódica y anódi- 
ca de cierre y apertura. Curvas intensidad-duración. 
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Abstract: This second part begins with the reformulation of bidomain theory of electrocardiology 
posed in the first part. The theoretical-analytical approach by means of nonlinear modal analysis is 
extended to take into account the unequal anisotropy of myocardium. Membrane recovery variables 
are taken into account employing a simplified mathematical model of excitation. The nonlinear 
dynamics of the electrode-myocardium system is studied just to threshold. The region of influence 
concept is re-examined and the generalized activating function is taking into account to construct 
the región of influence. A mixed homogeneous Dirichlet and nonhomogeneous Neumann boundary 
value problem is posed for the nonlinear partial differential equations of the bidomain model of 
electrocardiology. The polarizing effect of a pacemaker electrode is described, including virtual 
anodes and cathodes. Following the same general procedure implemented in the first part, the field 
equations are substituted by a system of nonlinear ordinary differential equations in the mode am- 
plitudes. In order to ease the calculation of the eigenfunctions and the eigenvalues of a Sturm- 
Liouville operator thar allows the definitions of the mode amplitudes, a virtual space (alibi) is 
employed such that in this space the ventricular myocardium can be represented approximately by a 
homogeneous and isotropic syncytium. The threshold value of a perturbation that causes instability 
is estimated by a modified version of Eckhaus method. A whole family of threshold space pat- 
terns of membrane polarization is identified in the state space constructed from the modal ampli- 
tudes, in the framework of the generalized mathematical model. The sizes and shapes of liminal 
regions (that must be depolarized over threshold of uniform membrane excitation in order to start a 
propagated action potential) are obtained for an anisotropic and heterogenous syneytium. Analytical 
expressions for cardiac strength-duration curves are derived, including closed form formulae for 
time constants (or chronaxies) and rheobases for cathodic make excitation, in terms of geometric 
and electrochemical electrode parameters, electrode-myocardium distance, electric properties of the 
volume conductor of the tissues, and the local distribution of electrotonic and excitability parame- 
ters of myocardium. Several experimental results, related with the unequal anisotropy of myocardi- 
um, are thus explained (liminal regions and the behaviour of the parameters of S-D curves) and 
some new facts are predicted (related with liminal regions, the family of electrode-myocardium 
threshold states and the parameters of the S-D curves). 


Key words: Bidomain models of the heart. Pacemaker electrodes design. Excitability thresholds. 
Unequal anisotropy. Heterogeneity. Alibi. Nonlinear modal analysis. Families of threshold states. 
Liminal regions. Cathodic and anodic make and brake stimulation. Strength-duration curves. 


(1)Introducción 


Como se dijo en la primera parte [20], el propósito principal de esta serie de trabajos es 
contribuir a una mejor comprensión de los problemas del umbral de estimulación eléctrica 
del miocardio, utilizando electrodos de marcapaso, a partir del análisis de un modelo 
matemático de la interacción entre el electrodo y el corazón. 

Los resultados de investigación aquí presentados desarrollan, amplían, corrigen y profun- 
dizan lo expuesto en la referencia [17]. El análisis se efectúa en el marco que brinda la 
teoría de sistemas [12] [13] [15]. 

En esta segunda parte, partiendo de la reformulación de los modelos de bidominio plan- 
teada en la primera parte, construimos un modelo matemático para analizar los fenómenos 
de umbral aplicable a un sincicio excitable no lineal, anisótropo y heterogéneo, con aniso- 
tropía desigual. La anisotropía desigual se refiere a la diferencia en el valor del cociente 
de conductividades paralela y perpendicular a la dirección local de las fibras miocárdicas: 
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en el dominio intracelular ese cociente vale aproximadamente 10 mientras que en el do- 
minio extracelular vale aproximadamente 2.5 [2] [3] [13]. Esta diferencia significativa 
permite comprender algunos de los mecanismos que subyacen a los fenómenos del um- 
bral de excitabilidad eléctrica del corazón, cuando la anisotropía desigual se considera 
junto con las características específicas de los canales iónicos en las membranas excita- 
bles del miocardio. 

Una de las consecuencias más notables de la anisotropía desigual es la diferencia en el 
mecanismo de excitación anódica de apertura del miocardio respecto del mecanismo de 
excitación anódica de apertura de las fibras nerviosas y de los tejidos excitables que se 
pueden describir mediante un modelo de sincicio eléctrico isótropo como el considerado 
en la referencia [20]. Debido a las propiedades de los canales de las células cardíacas, la 
excitación anódica de apertura que se observa en la práctica de los marcapasos parece 
depende de la migración de los cátodos virtuales que se generan debido a la anisotropía 
desigual, y no de los efectos debidos exclusivamente a las variables de recuperación de 
los canales [3] [5] [6] [21]. 

En esta segunda parte, una vez construido un marco general para el análisis modal, algu- 
nas de las características de los procesos de excitación vinculadas con la anisotropía del 
miocardio se deducirán y se discutirán aplicando herramientas de análisis modal no lineal, 
en forma análoga a la empleada en las referencias [18] y [19] para fibras nerviosas y en 
las referencias [17] y [20] para miocardio. Para tener en cuenta el efecto de las variables 
de recuperación de los canales iónicos, se empleará una modificación del modelo de 
Aliev y Panfilov [1] para la unidad de membrana de las células miocárdicas excitables. 
Asimismo, se considerarán, brevemente, algunas consecuencias de la heterogeneidad del 
tejido cardiaco. 


(2) Reformulación de los modelos de bidominio para el análisis del sistema electro- 
do-miocardio. 


Como se explicó en la referencia [20], para describir la polarización de las membranas 
excitables del miocardio debida a un campo de corrientes eléctricas producidas por un 
electrodo de marcapaso, es posible reformular las ecuaciones del modelo matemático de 
bidominio así: 


Lu (es Yalt.7)+ Jn Vu ») =Ve(Cevv,)+ve(DeJ) (3) 


0 

Ot 
Aquí 1 representa el instante de tiempo, mientras que F representa el vector de posición 
en el espacio euclidiano tridimensional. Los vectores se representan con una fleche sobre 
el símbolo correspondiente. La barra ondulada sobre un símbolo representa un tensor 
(díada) y V representa el operador gradiente. El punto grueso e representa un producto 
escalar de dos vectores, de un tensor (díada) con un vector, o de dos tensores. Como es de 


W =ElV,, 1) (1b) 
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uso común V .( ) representa el operador divergencia. El punto delgado - representa el 
producto de un escalar, un vector o una tensor por un número. 
La densidad de corriente eléctrica J (1,7) en el conductor de volumen de los tejidos, que 


aparece en el último término del segundo miembro de la ecuación (1 a), se supone es 
conocida. El continuo de las membranas excitables se describe mediante las variables de 


estado distribuidas (v,, (+,7),W(2,7)) siendo Vi (1,7) el campo de voltajes a través de 


las membranas excitables y W(t,7) el campo vectorial cuyas componentes son las varia- 


bles de activación y de recuperación de esas membranas. La evolución de dichas varia- 
bles viene dada por la ecuación (1 b). El número de variables involucradas (según el 


paradigma actual, asociadas a los canales iónicos) y la expresión concreta para F (Va, w) 


dependen del modelo de membrana utilizado, al igual que la expresión concreta para la 
densidad de corriente iónica a través de las membranas excitables, representada por 


JionlVrrW). En las referencias [2], [3] y [22] se pueden hallar diferentes modelos mate- 


máticos aplicables a la unidad de membrana de las células cardíacas. 
Con Cy se representa la capacidad por unidad de área de membrana y con z,y se repre- 


senta el área de membranas excitables por unidad de volumen de bidominio. 
Los tensores que aparecen en el primer y segundo término del miembro de la derecha de 
la ecuación (1 a) vienen dados por las fórmulas: 


C=f,-GjeGlef.-G, (Lo) B=elt..6,-1:6,) (1d) 


Los escalares f; y f. representan a las fracciones de volumen de los continuos intracelu- 
lar y etracelular (intersticial) respectivamente. Como se asume que el volumen ocupado 
por las membranas excitables es despreciable, resulta que f; +f. =1. Los tensores simé- 
tricos y definidos positivos G, y G, describen las propiedades de conductividad eléctrica 
de los continuos intracelular y extracelular, respectivamente. Por definición 
G= f. “G, +f.-* G, es un tensor definido positive y simétrico que caracteriza las propie- 
dades de conductividad eléctrica del tejido cardíaco (no obstante, ver la referencia [13] 
para una discusión de algunos aspectos delicados relacionados con el significado de este 
tensor) y G" es el correspondiente tensor inverso. De las formulas (1 c) y (1 d) se deduce 
que los tensores C and D son también simétricos y definidos positivos. 

La function activante generalizada Ve(Dej ) en el miembro de la derecha de (1 a) 


acopla el campo de densidad de corriente eléctrica en el conductor de volumen de los 
tejidos con las variables de estado del continuo de membranas excitables. 

En las escalas de tiempo involucradas en la operación de los marcapasos cardíacos los 
campos electromagnéticos son cuasi-estacionarios. Entonces el campo de densidad de 


corriente eléctrica en el conductor de volumen de los tejidos es solenoidal: VeJ=0 
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Como consecuencia, si los continuos intracelular y extracelular son homogéneos e isótro- 
pos, como se supuso en la primera parte de esta serie de trabajos, la función activante 
generalizada se anula idénticamente. El acoplamiento entre el campo de corrientes eléc- 
tricas producidas por el electrodo y las membranas excitables aparece a través de las con- 
diciones de frontera en la interfase que separa el miocardio de los otros tejidos [20]. 


Si los continuos intracelular y extracelular son isótropos pero heterogéneos, el tensor D 
se reduce al producto de una función escalar de la posición en el bidominio D(7) por el 


tensor unidad. Como consecuencia, la función activante se reduce a VD( r Jo J (r, 7). En 


ese caso, en un campo de corrientes no nulas, la función activante se anula solo si el vec- 
tor densidad de corriente es ortogonal al vector gradiente VD( r ) E 


Si los continuos intracelular y extracelular son homogéneos pero anisótropos, la función 
activante generalizada se anula idénticamente solamente para el caso de anisotropía igual. 
Pero si la anisotropía es desigual, como ocurre en el caso del miocardio ventricular, la 
función activante no se anula y el acoplamiento que resulta produce los denominados 
cátodos y ánodos virtuales [21] [22]. 

Como los campos de corrientes externos se asumen son cuasi-estacionarios, la densidad 
de corriente se puede expresar como el producto de un campo dependiente de la posición 


F ( r ) por la corriente eléctrica total ile Jinyectada a través del electrodo activo: 
J(t,7)=5(0)- E(7) 2) 
Entonces: Ve(DeJ)= ¡(1)-Ve(De F(7)= i(e).K(7) (3) 
El campo escalar K( r ) será denominado factor de forma de la función activante gene- 
ralizada o factor de forma de volumen: 
K(*)=V+(De F(7)) (4) 
Sea m7 ) el campo de vectores unitarios normales a cada punto de la interfase entre el 
miocardio y los tejidos externos. Suponiendo que en la frontera del bidominio la corrien- 
te eléctrica ingresa o egresa desde o hacia los tejidos circundantes solamente a través del 
continuo extracelular, en la superficie del músculo cardíaco la densidad de corriente en el 
continuo intracelular debe ser nula. En ese caso se puede demostrar que (ver [20]): 
ño(CoVVy )=ñelt, 6, e[G +3) (5) 
La ecuación (5) acopla, en la interfase entre el miocardio y los tejidos extra-miocárdicos, 
la polarización del continuo de membranas excitables con el campo de corrientes eléctri- 
cas producidas por el electrodo de marcapaso. Como no se anula idénticamente, aún si se 
supone que los continuos intracelular y extracelular son homogéneos e isótropos, el me- 
canismo representado por la condición en la frontera (5) es el único mecanismo de aco- 


plamiento entre el electrodo y las membranas excitables que siempre está presente. 
Para un campo cuasi-estacionario la ecuación (5) se reduce a: 


ñe(Covv, )=ñelr, O, e(G te F(5).i()=i(0)-S(—) 6 
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El campo escalar S (7), definido en la interfase entre el miocardio y los tejidos extra- 
miocárdicos, será denominado factor de forma de superficie: 

s(*)=5elf, 6, e(G* e F(+)) (7) 
Cuando el electrodo activo se encuentra localizado lo suficientemente próximo al miocar- 


dio (como es el caso de la estimulación a través de electrodos endocárdicos) la polariza- 
ción de las membranas excitables viene determinada fundamentalmente por la variación 


del campo F ( 7) en la ecuación (2) para la densidad de corriente en el conductor de vo- 


lumen de los tejidos. Esta variación actúa sobre las membranas excitables mediante dos 
mecanismos relacionados con el campo aplicado: uno es el que implica al factor de forma 
de volumen K( r ) definido en (4) y el otro es el que implica al factor de forma de superfi- 


cie S(7) definido en (7). 


(3) La región de influencia, las condiciones en su frontera y la identificación de un 
operador para el análisis modal no lineal 


Como se fundamentó en la primera parte, se necesita un dominio acotado B para poder 
desarrollar un análisis modal no lineal. Ese dominio acotado, denominado región de 
influencia del electrodo sobre el miocardio [20], se puede construir a partir de los facto- 
res de forma de superficie y de volumen. El factor de forma de superficie es una función 
escalar de la posición en la frontera del bidominio, mientras que el factor de forma de 
volumen es una función de la posición en el interior del bidominio. 

Considérense la familia de curvas (en la superficie frontera del bidominio) o la familia de 
superficies (en el interior del bidominio) en las cuales el correspondiente factor de forma 
(de superficie o de volumen, respectivamente) permanence constante. A continuación, 
considérense las líneas ortogonales a la familia de curvas o a la familia de superficies 
mencionadas, y tómese la restricción del correspondiente factor de forma, de superficie o 
de volumen, a una de estas líneas. Se obtienen así valores de los factores de forma K(r) (0) 


s(r) como funciones Ky(s) o Sy(s) de la longitud de arco s (abscisa curvilínea) corres- 


pondiente a la línea I' que se está considerando. Luego, se sigue la línea hasta abandonar 
la región donde el miocardio se encuentra significativamente polarizado. Se llegará a un 


último punto de inflexión S; de las funciones Ky-(s) oSy(s), luego del cual la restricción 


del correspondiente factor de forma tenderá monótonamente (y en general rápidamente) a 
cero. Calculamos la siguiente estimación de la distancia, a lo largo de la línea ortogonal 
TP, en la que la restricción del factor de forma sufre un decremento significativo hacia 
cero: 


Mm [Kp(s;)] 


[Sp-(s;)| 
s= 4 AA 
OO 


A (8 a) As= 4 (8 b) 
Sr) 
ds 
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Hecho esto, identificamos el punto de T' cuya abscisa curvilínea es s; + As y repetimos la 


operación para cada línea ortogonal posible. 

Cuando la función activante no se anula idénticamente, se define la frontera OB de la 
región de influencia como el conjunto formado por esos puntos de abscisas s; + As cons- 
truidos a partir del factor de forma de volumen, y la región de influencia B como el con- 
junto de puntos dentro de esta frontera. Así se dispone del dominio acotado que se bus- 
caba: lo bastante extenso como para involucrar los detalles espaciales asociados a la pola- 
rización de las membranas excitables, y lo bastante pequeño como para que la dinámica 
del umbral se pueda describir mediante un número relativamente pequeño de amplitudes 
de modo, como se hizo en las referencias [18] y [19] para el caso de la estimulación eléc- 
trica y magnética de fibras nerviosas. 

Cuando la función activante es idénticamente nula, se puede construir una región bidi- 
mensional, en la superficie del miocardio próxima al electrodo activo, formada por los 
puntos de abscisas s; + As determinados a partir del factor de forma de superficie, y se 
completa la región de influencia siguiendo el procedimiento explicado en la referencia 
[20]. Como en el presente artículo interesan las consecuencias sobre la dinámica del 
umbral de la anisotropía desigual y la heterogeneidad del miocardio, la región de influen- 
cia se construye a partir del factor de forma de volumen. Estas distinciones no se estable- 
cieron debidamente en la referencia [17], lo cual podría prestarse a confusión. 

La frontera OB de la región de influencia se puede dividir en dos superficies. Una de 


ellas, OB, , se ubica próxima al electrodo activo, en la interfase entre el miocardio y los 
tejidos extra-miocárdicos. La otra, OB, se ubica en el interior del miocardio, en la región 
de transición entre las membranas que están polarizadas y las que no lo están. 
Sobre OB), se aplica la siguiente condición de Neumann no homogénea: 
ño (Cev Vu )- i(e). s(r) (9) 

Si V, representa el valor de reposo del voltaje transmembrana y W, representa los valo- 
res de reposo de las variables de activación y recuperación de los canales iónicos, sobre 
OB, se aplica una condición de Dirichlet: 

Vu =Vx (10 a) Wu =Wz (10 b) 
Ahora, solamente falta identificar un operador lineal apropiado, cuyo dominio sea la re- 
gión de influencia, que se aplique a funciones que verifican las condiciones en la frontera 


(9) y (10), y cuyas funciones propias y valores propios se puedan utilizar en un análisis 
modal no lineal. Para identificar ese operador se re-escribe la ecuación (1 a) así: 


E —Yalt.F) == Sol We v el? 0 DY dio. xo( ) (11) 


Aquí R,, es una Resistencia de la unidad de área de membrana, 7y = Ry : Cy es una 
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A 52 Eg 5 . 
constante de tiempo de membrana, 2,,==%-.C es una constante de espacio tenso- 
XM 


a Ru As 
rial y Ko r )= va[ fu. pe F 6) es un factor de forma de volumen modificado. El 
XM 


factor de forma de superficie modificado es a su vez: Sol) A S (7) 

Xm 
Las ecuaciones (1 b) que describen la evolución de las variables de los canales iónicos 
permanecen incambiadas. 
La densidad de corriente iónica se desarrolla en serie de Taylor en torno a los valores cde 
reposo de las variables de estado de la membrana: 


Y OS; - OS; - is 
EV W)= ao We): (Vir =Ve)+ 2d (va W)e (im We)». 
Luego se introduce el siguiente operador lineal auto-adjunto: 


oJ ion % 7 
LVi]==Ru ¿7 (vo) (Uy; Vo) + Vo(2,,?e YY) (12) 
M 
Opera sobre funciones reales y regulares, definidas sobre la región de influencia B, y que 


verifican condiciones homogéneas de Dirichlet sobre OB, y condiciones de Neumann no 


homogéneas sobre 0B, . 


(4) Ecuaciones modales para la dinámica del umbral basadas en un modelo simplifi- 
cado de la unidad de membrana 


Si las variables de activación se consideran relajadas al equilibrio con el voltaje trans- 
membrana y las variables de recuperación se resumen en una sola, además del modelo 
clásico de Fitzhug-Nagumo (empleado fundamentalmente para fibras nerviosas) hay mo- 
delos específicos para describir efectos de umbral y propagación de potenciales de acción 
en el miocardio, como el de Aliev and Panfilov [1]. Estos autores discuten los inconve- 
nientes que presenta el modelo de Fitzhug-Nagumo cuando se intenta aplicarlo al miocar- 
dio y presentan un modelo matemático que permite superar la mayor parte de esos incon- 
venientes sin abandonar la descripción de la unidad de membrana mediante dos variables 
de estado. 

Para capturar los efectos asociados al umbral, es suficiente una aproximación polinómica 
a la densidad de corriente iónica al orden tres en v= Vy — Vr, como ésta, equivalente a la 


propuesta en la referencia [1], donde por definición w= W,yy — W es el apartamiento de 
la variable de recuperación respecto de su valor de reposo: 
Ru Sono, w)=v=b-v4c-v + a-v-w (13) 


Los parámetros b,c, «Y son positivos. En este caso la ecuación (11) se reduce a: 
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1on 


ca ole) —Rog + Sion lo, W)+V e (4,2 Vo)+ i(e)-K,(7) (14 a) 


Para los fines del presente trabajo se asumirá que la cinética de la variable de recupera- 
ción viene dada por la fórmula siguiente, que es una modificación de la correspondiente 
ecuación del modelo de Aliev y Panfilov: 


1) pS C 2 
T,:—WMt,r)=6o-—-v -W 14b 
"a (t,7) a Y (14 b) 


El parámetro y es positivo y 7,, es una constante de tiempo que por lo general es un orden 
de magnitud mayor que la constante de tiempo de membrana 7 yy . 
El operador de Sturm-Liouville (12) se puede escribir ahora así: 
Lo]=Ve(42eVo)+v (15) 
Sus funciones propias, definidas en la región de influencia del electrodo y verificando 
condiciones homogéneas mixtas de Dirichlet y de Neumann en su frontera, son ortogona- 
les dos a dos y se pueden hallar resolviendo la ecuación de valores y funciones propias: 
7 2 2 

vela, *Vo,)=4, 0, peliZoms (16) 

Siguiendo el procedimiento descripto en las referencias [14], [17] y [20], empleando 


funciones propias normalizadas, introduciendo los parámetros A, =1+ xo y las solucio- 

nes tentativas VL=>A, ()-p, (7) y W= 2B, ()-p, (7), se obtienen las ecuaciones para 
P q 

las amplitudes de modo, (con p=1,2,...): 


dA E E 
y 
Tu * =-A, Ay +b* A A ge A ASE 
dt q.r=1 q.r,s=1 


o “As B, + f, (1) 
qg.r=1 


(17 a) 
dB C E 
B.= 
crab: ra Lips Ay B=1 Br (17 b) 
q.r= 
En estas fórmulas, por definición: 


Doar = Pp" Pq" Pr: AV (18 a) Irars = Op" Py O," Ps: AV (18 b) 
B B 


T 


Teniendo en cuenta las expresiones para los factores de forma modificados 


a Ru ALS ZA R Lo ; , 
Ko(7 Java Lu. 5er) y Sol7)== . s(7), los términos forzantes vienen dados 
XM XM 

por las fórmulas (para p=1,2,3,...): 
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FO = 40) [op rav + Jo, (7) sol) as (19) 


OB, 
El término forzante de volumen f,, y (1)= 10 Jo, (o) ar] es la proyección del 
B 


factor de forma de volumen sobre la función propia normalizada a (7 ) y permite estu- 


diar los efectos de los cátodos y ánodos virtuales sobre la dinámica del umbral. 


El término forzante de superficie f, (1) = ile). lo a (7). So (+). ds es la proyección del 
OB, 


factor de forma de superficie sobre la función propia normalizada p,(r) , y está siempre 


presente y operando para polarizar el tejido durante un pulso de corriente eléctrica aplica- 
do, aún si el sincicio se considera homogéneo e isótropo. 

En la primera parte de este trabajo se mostró que para un sincicio homogéneo e isótropo 
es posible estudiar la dinámica del umbral para pulsos catódicos (excitación catódica de 
cierre) empleando una amplitud de modo dominante (amplitud de primer orden A, ) co- 


rrespondiente al voltaje trans-membrana. Para estudiar la excitación anódica de apertura 
en un sincicio semejante y con un mínimo de complejidad, se requiere dos amplitudes de 
modo dominantes: una correspondiente al voltaje trans-membrana A; y la otra correspon- 


diente a la variable de recuperación B;,. 


Cuando el sincicio presenta una anisotropía desigual muy marcada, como por lo general 
ocurre con el miocardio ventricular, el escenario es más complicado debido a los efectos 
de los electrodos virtuales [21]. 

Esto debe ser tenido en cuenta al desarrollar un análisis modal no lineal, trabajando con 
amplitudes de modo de orden superior al primero y empleando al menos una variable de 
recuperación en el modelo matemático de la unidad de membrana. En todos los casos el 
conjunto de amplitudes de modo utilizadas es finito: se desprecian las de orden superior al 
máximo necesario para capturar los aspectos interesantes de la dinámica. 


Reachable 
Sets from 
the origin 


Pseudo-Threshold manifold 


Figura 1. Esquema de una partición del espacio de las amplitudes de modo. 


Rev. SCP” Volumen 19 N*2 2014 pp. 135-168 ISSN 0379-9123 


Los resultados de diferentes casos investigados por análisis modal no lineal, tanto de 
estimulación eléctrica o magnética de fibras nerviosas o musculares esqueléticas [14] [18] 
[19], como de estimulación eléctrica del corazón [14] [20], presentan ciertas característi- 
cas comunes que se pueden resumir en el croquis que aparece en la Figura 1 precedente, 
tomado de la referencia [17]. 

Para cada sistema electrodo-tejidos-elemento blanco de la estimulación, el espacio de las 
amplitudes de modo se puede dividir aproximadamente en dos subconjuntos: el conjunto 
de decaimiento y el conjunto amplificador. El estado de reposo de las membranas, repre- 
sentado por el origen del espacio de las amplitudes de modo, pertenece al conjunto de 
decaimiento. También pertenecen a este conjunto aquellos puntos tales que si el estado 
del sistema se ubica en ellos al final de un pulso estimulador, entonces el sistema retorna- 
rá al estado de reposo sin presentar un potencial de acción. 

Si el pulso deja el estado del sistema en el conjunto amplificador, emerge un potencial de 
acción y la estimulación es exitosa. 

Los estados fronterizos entre los puntos del conjunto de decaimiento y los puntos del 
conjunto amplificador se pueden considerar como estados umbral: para cada sistema 
electrodo-tejidos-elemento blanco hay una familia de estados umbral. 

En los modelos cuya unidad de membrana involucra variables de recuperación, la frontera 
entre el conjunto de decaimiento y el conjunto amplificador es más o menos difusa: en 
este caso hay una pseudo-variedad umbral, indicada en la Figura 1. 

Además, esta figura esboza los subconjuntos del espacio de estados alcanzables desde el 
origen (desde el estado de reposo) mediante la aplicación de pulsos catódicos o anódicos. 
Se observa que ambos conjuntos alcanzables parten del estado de reposo, atraviesan el 
conjunto de decaimiento, cortan la variedad umbral y se internan en el conjunto amplifi- 
cador. 


(5) Alibi, factores de forma, electrodos virtuales, regiones de influencia, funciones y 
valores propios para un modelo del miocardio ventricular 


Para hallar las regiones de influencia del electrodo de marcapaso sobre el miocardio, para 
determinar las funciones y los valores propios que se utilizan para deducir las ecuaciones 
de evolución de las amplitudes de modo, y para poder calcular las proyecciones de los 
factores de forma sobre las funciones propias, es necesario tener en cuenta la anisotropía 
desigual y por lo menos cierto grado de heterogeneidad en el tejido. 

Todo esto se puede abordar por métodos numéricos ab-initio, empleando algunos de los 
poderosos paquetes de software comercial disponibles. No obstante, puesto que el propó- 
sito de esta serie de trabajos es utilizar métodos analíticos (hasta donde esto resulte facti- 
ble), en lo que sigue se dará una solución analítica aproximada al problema de las funcio- 
nes y valores propios y al problema de determinar los factores de forma. 

Teniendo en cuenta las dimensiones que poseen los electrodos de marcapaso actualmente 
utilizados para la estimulación del miocardio ventricular, las curvaturas de la superficie 
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del endocardio y las características de anisotropía y heterogeneidad de este músculo, es 
posible introducir un espacio virtual (alibi) en el cual un problema tridimensional comple- 
jo de polarización de las membranas excitables por un electrodo de marcapaso activo se 
puede aproximar por un problema más simple de polarización en un medio isótropo y 
homogéneo, resoluble por métodos analíticos. Como entre los puntos del espacio virtual y 
los puntos del espacio físico hay una correspondencia regular uno a uno, una vez resuelto 
el problema en el espacio virtual, se puede retornar a las coordenadas originales en el 
espacio físico donde se pueden evaluar las modificaciones de la polarización de las mem- 
branas excitables debidas a la anisotropía y la heterogeneidad. 

Además, el método del alibi resulta aplicable a la investigación tanto de la propagación 
del potencial de acción [14] como de la sensibilidad de detección de señales eléctricas por 
los electrodos de marcapaso [16]. 

Se supone que un pequeño electrodo activo se halla ubicado en la cavidad ventricular, 
tocando el endocardio. Para los fines del presente análisis la superficie endocárdica se 
puede representar, localmente y con suficiente aproximación, mediante un plano. Las 
fibras miocárdicas corren, casi paralelas, en capas superpuestas. La dirección de las fibras 
de cada capa varía regularmente desde el endocardio al epicardio. 


Se introduce un triedro Ty de vectores unitarios (versores) doc Toda con el origen y los 
vectores iy, jp fijos en el plano endocárdico. El versor ko es orthogonal a ese plano y se 
elige apuntando hacia el epicardio. Se considera la semirrecta definida por el versor ko y 
el origen del triedro T¿. Cada punto de esa semirrecta, de coordenada z , se toma como 
origen de un nuevo triedro ortogonal T, de versores €, =€y, €), €, = ko , donde e, da la 


dirección de las fibras en la capa situada a una profundidad Z a partir del epicardio. En- 
tonces, cada punto en este modelo geométrico de miocardio ventricular queda determina- 


do por una coordenada z y por sus tres coordenadas X, , X, , X, respecto del triedro local 
T, . Entre las coordenadas X, y relativas a TI, y las coordenadas x,x, relativas a T, se 
verifica la relación: Xx =x: sinlo(z)]+ y: cos|a(z)] Q0 a) 
X) =-X* cos[ar(z)]+ y: sin[a(z)] (0 b) 
En estas ecuaciones alz) es el ángulo comprendido entre el versor e; (que da la dirección 


de las fibras en la capa considerada) y el versor do a 


(5.1) Cálculo del factor de forma mediante el alibi para un modelo de miocardio con 
anisotropía desigual pero heterogeneidad débil. 


Si G es el tensor de conductividad del conductor de volumen equivalente a los tejidos y 
D es el correspondiente potencial eléctrico, la densidad de corriente viene dada por la 


relación J = -G eVO Como la divergencia de J es nula, el potencial verifica: 
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ve(Gevo)=0 21) 
El tensor de conductividad es diagonal respecto de T, en cada capa de fibras miocárdicas, 


pero sus direcciones principales varían con continuidad de una capa a otra, como se puede 
deducir de las ecuaciones (20 a) y (0 b). 
Supóngase que, con suficiente aproximación sus componentes diagonales dependen so- 
lamente de la coordenada principal correspondiente, es decir vienen dadas por funciones 
del tipo Gi (x,), Gr (x> ),G3(x3). Entonces, si D(t, x,,x,,x3) es una representación del po- 
tencial eléctrico en las coordenadas locales asociadas a T_ la ecuación (21) se reduce a: 
3.0 0 
YN G,(x,)—0 |=0 (2) 
n=10X,, X, 
No obstante, todavía permanecemos en el espacio físico, rotulando a los mismos puntos 
de otra manera (alias, mismo punto, diferentes coordenadas). Ahora se introduce un espa- 
cio virtual (alibi) cuyos puntos vienen caracterizados por las coordenadas Ele da 5s , que 
se encuentra ocupado por un conductor de volumen homogéneo e isótropo de conductivi- 
dad G, . Definimos la conexión entre los puntos correspondientes de uno y otro espacio 
mediante las fórmulas: 


Xx Go 
E=S a e=123 (23) 
0 ¿Nu 


Gilu) 
La conexión entre la representación local del potencial en el espacio físico D(t, x,,x,,x3) 
y el potencial en el punto correspondiente del espacio virtual wlt,E,, Es, Es) se define por 
la relación: 
Dlr,x;,x2,3)=y (1, E (1), 6, (1) 5(13)) (Q4a) 
Entonces, la corriente en el espacio virtual i,(1) y la corriente en el espacio físico ie ) 


a el 
deben estar relacionadas así”: 


40 ES 0) As 


Si el potencial en el espacio físico D(t, x,,x,,x3) verifica la ecuación (22), entonces el 
potencial wle, E, Ez) en el espacio virtual verifica la ecuación: 
20% E 
yy PA (25) 
k=1 05% 2 Lock 


En la ecuación (25) aparecen dos escalas espaciales en el espacio virtual. Una de 


! Una justificación de esta relación puede hallarse en la solución de un problema electroestático 
análogo, planteado en un medio anisótropo y homogéneo. Ver L. Landau y E. Lifchitz, “Electrody- 
namics of continuous media”, Pergamon Press, N. Y., 1960, pág. 61. 
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Oye 


0 
ellas, L,, y; % E 
OYw 
06; 
el gradiente del potencial y sufre una variación significativa. 

oa du G; (x, (E, )) 

ROS 
| dG; x (Es )) 
ds, 


través de la cual la componente principal G, del tensor de conductividad sufre una va- 


estima una distancia a lo largo de la coordenada ¿, a través de la cual 


La otra, estima una distancia a lo largo de la coordenada Ce a 


riación significativa. Si L£,,, es al menos un orden de magnitud inferior a Lg, (es decir, 

si la heterogeneidad es débil) el potencial en el espacio auxiliar wlt, EE, E) verifica 
aproximadamente la ecuación de Laplace: 

3 yw 7 

OE? 


Para un semi-espacio con una corriente total i, (1) inyectada en el origen en condiciones 


(26) 


cuasi-estacionarias, la ecuación (26) posee la solución: 
1 


AO) 
wle, Ei Es 67) 27 Go : ler ere 


Entonces, de (24 b) se desprende que el potencial en el espacio físico, cuando la hetero- 
geneidad es débil, viene dado por: 


¡ 1 
DÁt, x Xx )= ile) E (28) 
EAS 277 -,[G,(x, )-G,(x,)-Gy(x3) ¡ xi x x 
G; (+) Ga(x,) Gs(x3) 
Cuando el miocardio se puede considerar homogéneo, la expresión (28) es exacta. Si, 


(27) 


además, el hemi-espacio considerado es isótropo, de modo que G,=G,=G;¿=G, la 


fórmula (28) se simplifica, como debe ser, a la expresión bien conocida: 


LO) 1 


Para simplificar, para el cálculo del factor de forma se supondrá un continuo anisótropo 
pero homogéneo. En ese caso el factor de forma de volumen (ecuación (4)) se reduce a 


OF, (1,7,2) 
OXy 


0 J¿lt,x1,x2,x3) 


2 
e 3 Md 


3 3 
K(x,,x2,x Je YD, 4 = YD, é 
k=1 k=1 
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A partir de (1 d), (2) y (29), teniendo en cuenta que J, =-G; —0 se obtiene: 
Xx 

3 0 dr x,x,x 

2D GE 2 y La 3) 

k=1 OXz Ni 


OF AX, 2,3) _ 
OXz 


3 
K(x¡,x7,%3 )= 24 


DG, == (fo: —£, is) 


A partir de estas últimas fórmulas y de (28) escrita para un continuo homogéneo, resulta, 
después de varios pasos de cálculo: 


ls x 0 Ss na E 
dis Eg 3" 
K(x X2>X3 Je a ; % S Sa (20) 
es EC Es 
CE (1, 2,3P 
G G, G; 


Por definición: 
po =D,+D3-2:-D,  ,=D,+D;,-=2:D, 3 =D, +D,=2-D, 


Entonces (+, +43 =0, de modo que cuando la anisotropía es desigual, en general 


hay siempre al menos un gp positivo y al menos un gy negativo. Como en el caso del 
miocardio ventricular se pueden efectuar las aproximaciones G,,=G,3 y G¡2=G;3 


([31,122]) resulta G, =G;z y D, =D, con lo cual q, =2-(D, =D) y P=0=D,-D,. 


Entonces, con p= D, — D,, el factor de forma de volumen se reduce a: 


G 
o [6 (A 
K(%,0,%)= 2 | 7 y 
1 


5 
2, 21,6 21 
FJ A 
( 2 2) E 1 j] 
(5.2) Construcción de la región de influencia del electrodo y caracterización de los 
electrodos virtuales. 


(31D) 


Se puede construir una región de influencia de un electrodo de marcapaso sobre el mio- 
cardio, aplicando el procedimiento sugerido en la parte 3, ecuaciones (8), a la fórmula 
(31) para el factor de forma, suponiendo que el electrodo real se sustituye por un electro- 
do puntual ubicado en su centro eléctrico. La posibilidad de sustituir un electrodo de mar- 
capaso por un sistema equivalente de multipolos de corriente ubicados en el centro eléc- 
trico se estudia en las referencias [13], [14] y [16]. 

Como se muestra a continuación, la región de influencia, construida a partir del factor de 
forma de volumen, resulta ser un hemi-elipsoide de semiejes a; (paralelo a las fibras en el 


plano del endocardio), a, (ortogonal a las fibras en el plano del endocardio) y az (orto- 
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gonal al plano del endocardio y en dirección al epicardio), situado a una distancia d del 
punto donde se ubica el centro eléctrico del electrodo: 


2 2 2 
X; x x-d 

+ EZ na (32) 
4 4 4 


Las coordenadas x, están referidas al triedro local T_ que va cambiando de orientación al 


pasar de una capa de fibras a otra. 
Sustituyendo las coordenadas respecto a T, por las coordenadas x, y,z respecto del 


triedro fijo triedro T, , mediante las fórmulas (20 a) y (0 b), se obtiene como región de 
influencia un hemi-elipsoide sometido a una torsión variable en torno al eje z=x3 dada 
por el ángulo alz). 

Si un electrodo puntual se encuentra situado en el origen y el tejido excitable forma un 
hemi-espacio ubicado en xz > d , la restricción del factor de forma de volumen a la super- 


(33) 


ficie endocárdica es: K(x,,x2,d )= P [E . 


La coordenada Xx, es paralela a las fibras, mientras que x, es perpendicular a las fibras. 


e 
em eS Qp G G; 
Tomando x, =0 resulta: K(x,,0, d) 2. (34) 
Gi 


(2. +4) 
Gi 


27 
Es equivalente a la función activante unidimensional, para una fibra nerviosa o muscular, 
cuya gráfica es semejante al perfil de un sombrero mexicano. (En el marco de la neuro- 
estimulación por electrodos externos fue introducida por Rattay [21]). 
Dependiendo del signo del pulso de corriente, K(x, 0,d ) representa en el intervalo 


G G z % j 
—L ds x £ —L_ .d ya sea un cátodo, ya sea un ánodo, mientras que a uno y 
2-G, 2-G, 


N]|u 


otro lado de este intervalo, a lo largo de la abscisa X, aparecen ya sea un par de ánodos, 
ya sea un par de cátodos virtuales, respectivamente, antes de que K(x, 0,d ) se aproxime 


monótonamente a cero al tender x, a +00. 


Si x,,es el punto de inflexión positivo más alejado de x, =0, definimos así la longitud 
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K(s,,0,4] fa, 


del semieje paralelo a las fibras: A=X ¡+ = -d 35 
Ip Fx 0,4) 16, e 
Ox, 
Tomando x, =0 se tiene: K(0, x,,d )= e: 16 : : E (36) 
27 V G; d 2) 
xx +d 


Dependiendo también del signo del pulso de corriente, K(0, x2,0 ) ya sea representa un 
cátodo, ya sea un ánodo, a lo largo de la abscisa x,. En este caso no aparecen electrodos 
virtuales. Si x) es el punto de inflexión positivo de (36), definimos así la longitud del 


semieje ortogonal a las fibras y contenido en el plano del endocardio: 


Ñ [K(0,1,1,4)] 
4) = X2 ¡+ Ko. a =d 87) 
Ox» 


La Figura 2 muestra un croquis conjunto de K(x, 0,d ) (restricción del factor de forma 
sobre el eje paralelo a las fibras) y K(0, Xx7,0 ) restricción del factor de forma sobre el eje 


perpendicular a las fibras). 


ss Fl%z, A dy 


Flo, ER 


Figura 2 


Además, figura 2 muestra la determinación de los semiejes a; (paralelo a las fibras en el 
plano del endocardio) y a, (ortogonal a las fibras en el plano del endocardio). 


Teniendo en cuenta ambas abscisas a un tiempo, se obtiene el patrón de líneas de nivel del 
factor de forma que se esboza en la Figura 3. Obsérvese la deformación de las líneas de 
nivel que caracterizan al denominado cátodo o ánodo en forma de “hueso de perro” en el 
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plano del endocardio. A los costados, sobre el eje paralelo a las fibras se dibujaron los 
contra-electrodos virtuales (ánodos o cátodos). La separación entre líneas de nivel se ha 
dibujado solamente para fines ilustrativos y no representa una misma diferencia de poten- 
cial. La débil asimetría que muestra el croquis sugiere lo que podría ocurrir debido a una 
heterogeneidad débil en el miocardio ventricular. 


Figura 3 


En la referencia [3] se pueden encontrar resultados sobre el patrón de polarización “hueso 
de perro”. Algunos fueron obtenidos mediante experimentos con corazones en modelos 
animales, utilizando trazadores fluorescentes para los cuales la permeabilidad de las 
membranas excitables varía con el grado de la polarización eléctrica. Otros resultados 
fueron obtenidos a partir de simulaciones digitales aplicando métodos numéricos ab-initio 
a partir de modelos matemáticos discretos que representan un tejido perfectamente homo- 
géneos pero con anisotropía desigual. 

La restricción del factor de forma al eje perpendicular al plano del epicardio es: 

p [6 1 


K(0.0,1,)=>= e (38) 
1 3 


En este caso la restricción del factor de forma disminuye monótonamente al aumentar la 


coordenada x, > dl . El semieje az del hemi-elipsoide de influencia se puede estimar así: 


[K(0,0,d)| _ 4 
ad =2.d (39) 
óK(0.0,d) 3 
OXz 


(5.3) Determinación de las funciones y valores propios mediante el alibi para un 
modelo de miocardio con anisotropía desigual pero heterogeneidad débil. 


Tomando en cuenta que la constante de espacio tensorial se relaciona con el tensor 
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C= f, -G, .G7" ef, -G, por medio de la ecuación iS = En 6 , la ecuación (16) se 
XM 
Td E Ry = 2 
puede re-escribir así: Ve -CoVo, |=-4, Op (40) 
XM 


Debido a la simetría de los tensores de conductividad intracelular G, y extracelular G, 
(121 [5] [6] [21]) se puede asumir que el tensor C =f; :G .G7 ef. «6, es diagonal res- 
pecto de T, en la capa correspondiente de fibras miocárdicas. Con suficiente aproximación 
sus componentes principales C;,C,,Czse pueden representar como funciones de una sola 
coordenada C; (és) ,C) (x,) y C; (x,). Entonces, si Y, (2128) es la función propia en 
coordenadas locales, (40) se reduce a: 

3 

¿qu :C, he) =—4p “Pp (41) 
Ahora, introduzcamos un nuevo espacio virtual (alibi) de coordenadas 7),27,,773 , lleno 


con un sincicio eléctrico homogéneo e isótropo de constante de espacio escalar 


[R = pa 
Am =. LC, , con un escalar C, en lugar del tensor C. La conexión entre puntos 
Xm 


correspondientes de uno y otro espacio viene dada por las ecuaciones: 


Xx Co 
My = - ds k =1,2,3 (42) 
oy (s ) 
Definimos ahora Yin» 721)=09, Lar). 2(72).307,)) (43) 
Introducimos las funciones: Ch 77 ) =C; lx, (77 ) k =1,2,3 (44) 
A partir de (42), (43) y (44) se desprende, luego de varios pasos de cálculo, que (41) se 
302 
puede re-escribir así: A y Yo 1+ z : 2 —u5 Y) (45) 
k=1 ON; 2 Lor 
0 
e LL (46 a) L, y (46 b) 
y, ? AC 
One Um 


Si las escalas L,, , son al menos un orden de magnitud menor que las escalas £, , (hete- 


rogeneidad débil), entonces las funciones propias y y (17,172,273) en el espacio auxiliar 
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2? 
verifican aproximadamente: A y E = 15 Y, (47) 
k=1 ON; 


L, 
La aproximación mejora cuando ec disminuye. 
c,k 
A la región de influencia B del electrodo de marcapaso sobre el miocardio en el espacio 
físico, le corresponde una región B, en el espacio virtual. A las fronteras 0B, (donde se 
anula la derivada en dirección normal de la función propia ,) y 0B, (donde se anula 
la función propia Q,) les corresponden las fronteras 0B, , (donde se anula la derivada 
en dirección normal de la función propia y, ) y 0B), (donde se anula la función propia 
Y ,). Considérese ahora un sincicio con anisotropía desigual, pero homogéneo. A la 


región de influencia definida por las desigualdades (32) (con respecto de las coordenadas 
del triedro local T, ), le corresponde en el espacio virtual la región de influencia definida 
por las desgualdades: 
2 2 2 
15 2 y de) Ly mzd, (48) 
bi bz b; 


La región de influencia en el espacio virtual es también un hemi-elipsoide de semiejes: 


b, = [Sa k=1,2,3 (49) 
k 


El plano que corresponde, en el espacio virtual al endocardio, se halla a una distancia del 


centro eléctrico del electrodo dada por: d, = e -d (50) 
3 
De las estimaciones (35) para a; , (37) para a, y (39) para az, se obtiene: 


bz [Soy [Sd (Sl a) b, = ¡ey (51 b) os ed (31 c) 


El problema de valores y funciones propias (47) con condiciones homogéneas mixtas en 


UN e ¿ld -1 
2 


la frontera del hemi-elipsoide (48) (de Dirichlet sobre z 7 con 
bi bz b; 
2 2 
1), 2 d,, y de Neumann sobre mm + pa =1 con 7¿=d, ) se puede resolver en términos 
1 2 


de funciones de Lamé separando variables mediante coordenadas elípticas [8] [11]. De- 
bido a la simetría espacial que existe en este caso, las soluciones que buscamos son las 
mismas que se obtienen resolviendo el problema sobre el elipsoide completo 
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2 a 2 
=d Ps , em 
de +B + Cr > .) <1 con condiciones homogéneas de Dirichlet sobre su frontera 
bi bz b; 


mm lm-d.) race 
DAR a =1. No obstante, para los fines de la presente investigación, es 

bi bz b; 

suficiente aplicar los resultados presentados en la referencia [20] para acotar los valores y 


las funciones propias de un problema del tipo vu, =v,? -u, en un dominio (2, con 
condiciones homogéneas de Dirichlet sobre 0(2. Es posible acotar inferior y superior- 
mente los valores propios recurriendo a regiones para las cuales el problema posee una 
solución analítica conocida. Supóngase una familia de regiones Q(s), donde £ es un 
parámetro real (o un conjunto de tales parámetros). Resolvemos el problema 
—V?u(F)=v? -u(r) sobre Q(s) con la restricción u(7)=0 para todo F coQls). Se 
obtiene un conjunto de funciones propias u,, (F; 2) y de valores propios correspondientes 


2 s á X > y 
v (s). Si las dimensiones de las regiones dependen en forma regular del parámetro e, 


lo mismo ocurre con los valores propios v. (€) y con las funciones propias u,(7;8). 
Entonces es posible demostrar (ver, por ejemplo, la referencia [11]) que si 
O(s,)= O(s)= O(s)), donde la inclusión es en sentido estricto, los valores propios co- 


rrespondientes al mismo índice n (o índice múltiple /,m,n) verifican la desigualdad 
estricta v,?(s,)>v, (€) >v, (£2)- Por ejemplo, si Q(£) es un elipsoide de semiejes a, 


by c podemos tomar como región Q(s,) el paralelepípedo rectangular con el mismo 


2a 2b 2:c. sa a ) 1 (étoed ; 1 
Aa y a y como región Ols,) el paralelepípedo rectangular 
análogo pero con lados 2-a, 2-by 2-c. Como puede verse en las referencias [8] y [11], 
si (2 es un paralelepípedo rectangular de lados 2-a, 2-b y 2-c los valores propios son: 


2 2 2 2 
2 T l m n 
Viomn ==" + 52 
l,m, 4 [5 p? .) ( ) 


centro y lados 


Las funciones propias son: 


Lx m-7:x Nn*T-X 
,y,z7)= sI - si - si 53 
ala y,2) al a ) sl 7 ) al 7 ) (53) 


Entonces el valor propio Vimiale) del problema planteado sobre el elipsoide verifica 


dE) < Yoma (e)< Vandal) , es decir: 


2 2 2 2 2 2 2 2 
TT l mon 2 3-7 l mon 
( 2 + 2 emo (s)< ( + + - (54) 
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(6) Regiones liminales y curvas intensidad-duración cuando el primer modo se puede 
considerar dominante 


Como una primera aproximación a los problemas que plantea el umbral de estimulación 
del miocardio, teniendo en cuenta la anisotropía desigual, se asumirá que el primer modo 
domina la dinámica. Como consecuencia es posible desacoplarlo de los modos de orden 
superior, y entonces las ecuaciones (17 a) y (17 b) se reducen a las siguientes: 


A; 
TO E (55 a) 
dB; C 2 
Ty" =A, -A -B 55 b 
ds E YB (55 b) 
hu=50-0-p:dV  (55c) hiua=50 0:00 :4V (55d) 
B B 


Las integrales 1,¡¡y f,¡,, Son positivas porque q, , que no cambia de signo en el interior 


de la región de influencia B del electrodo sobre el tejido excitable, se toma positiva. 
El parámetro geométrico K, y el coeficiente sin dimensiones 4, vienen dados por: 


K,=|pj(7) Ko(F).aV+ [o (7)-Sy(F)- ds. (55e) y A=l+m (551) 
B OB, 

En (55 6, y i es el valor propio positivo correspondiente a la función propia normaliza- 

da de primer orden (, del operador de Sturm-Liouville introducido en la parte 4 de este 


trabajo. 
Los factores de forma modificados que aparecen en (55 c) se relacionan con los factores 
de forma de volumen y de superficie mediante las ecuaciones ya vistas en la parte 4: 


Ko(7)=v9( 22-57) ys re sr) 


M 


Cuando se anula el pulso de corriente aplicado, el único punto de equilibrio de (55 a) y 
(55 b) es el origen (el estado de reposo de la membrana excitable). Un análisis análogo al 
efectuado en la referencia [19] para el modelo de membrana de Fitzhugh-Nagumo, permi- 
te concluir que el plano de fases se puede dividir en un conjunto de decaimiento y un 
conjunto amplificador, separados parcialmente por una curva umbral. Esa curva se puede 
aproximar por una barrera umbral recta, como la que se sugiere el croquis de la Figura 4. 
Los estados en principio accesibles por pulsos de corriente catódicos y anódicos, aplica- 
dos a través de electrodos de marcapaso, se indican como regiones grises en el plano de 
estados. Tanto el conjunto de estados accesibles mediante pulsos catódicos como el con- 
junto de estados accesibles mediante pulsos anódicos atraviesan el conjunto de decai- 
miento a partir del origen y alcanzan el conjunto amplificador. 

(El análisis detallado de la dinámica, análogo al efectuado para una fibra nerviosa en la 
referencia [19] se presentará en otro trabajo). 
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brad 


barrera um 


EAN cos 


Figura 4 


(6.1) Regiones liminales para el modelo de miocardio con anisotropía desigual. 


Como puede verse en el croquis de la Figura 4, la barrera umbral es tangente, en uno de 
los dos puntos donde esa curva corta el eje A, , a la curva: 


1 2 
B, = A A (56) 
a ti11 


De (36), se obtiene para la barrera umbral la fórmula: —B,= amo.) . (A, ES A.) (57) 
1 
Ese punto Aj, , si existe, es uno de los puntos de equilibrio de la ecuación que describe la 


dinámica de la membrana cuando la variable de recuperación permanece fija en su valor 
de reposo B,=0 y el pulso estimulador ya ha cesado: 


dA, 
a TACA bl AR Ci AS (58) 
El estado de dE A,=0 siempre es un punto de equilibrio estable. Si el discriminante 
BL E! (1 +u 2 ) es positivo (tener en cuenta (55 f) para relacionar A, con 
1 
META CL 

e ), hay otros dos puntos de equilibrio, uno inestable y — bl JA Y Otro esta- 

al 111 

bl ; a 
ble A, == ——+ A - (Esto, y sus consecuencias sobre la dinámica del umbral, ya fue 

E 2c155 


estudiado en detalle en la referencia [20] para un sincicio eléctrico homogéneo e isótro- 
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po.) 


Figura 5. Diagrama cualitativo que muestra una bifurcación global de las amplitudes de equi- 
librio del modo dominante en función del valor propio del modo dominante. 


La Figura 5, tomada de la referencia [20], muestra un diagrama de bifurcación de las tres 


amplitudes de equilibrio en función de un parámetro proporcional a 14, . Reinterpretando 


el significado de 4, , el diagrama resulta aplicable cuando el sincicio funcional es anisó- 


tropo y heterogéneo. 
La línea punteada vertical en la Figura 5 asciende desde un valor crítico del parámetro P, 
2 2 
ESTAN 


que corresponde al valor crítico? Mic = 3 1 del parámetro yu o , hasta el 
É 4-c 


punto donde se produce una bifurcación global del umbral. Para ese valor de y io se anu- 


la el discriminante A. Para valores mayores del parámetro el umbral desaparece. Para 
valores menores, además de la rama coincidente con el eje de abscisas y correspondiente 
al estado de reposo, se tienen dos ramas adicionales: la inferior corresponde al umbral y la 
superior a un estado excitado estable que resulta posible cuando no hay recuperación de la 
membrana de las células excitables. 

Cuando el tejido presenta anisotropía desigual y puede considerarse en primera aproxima- 
ción homogéneo, para las regiones de influencia y los valores propios analizados en las 
partes 5.2 y 5.3 de este trabajo, asumiendo iguales propiedades eléctricas en las direccio- 
nes ortogonales a la dirección de la fibra tanto en el continuo intracelular como en el 


j , R , , 
continuo extracelular, el valor propio Mo= —.C, ei del modo dominante, aplican- 
M 


? En lo que sigue se supone que ese valor crítico es siempre positivo. 
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R 
do la ecuación (54) queda acotado inferiormente por 2..c,.—-. + + y 
XmM 4 b; 


E 3771 1-1 
superiormente por «Co: AGS 
Xm 4 bE by Dd, 
Con suficiente aproximación para los propósitos del presente trabajo y teniendo en cuenta 


las estimaciones (51 a), (51 b) y (51 c) de las dimensiones b,, b, y bz del hemi- 


elipsoide de influencia obtenemos finalmente: 


Se G, 5 1 
Mz o A LN (59) 
8 Xm G 4 d 
Entonces la condición para que el sistema sea excitable puede resumirse así siendo 
E R G, 5 1 Buf 
sa 1 - 111 di Mola PA El: sE de (60) 
á 4.c Am G 4 d 5-7 


Como se hizo ya en la referencia [20], se denominará región crítica o región liminal a 
una región tal que una despolarización en una región estrictamente incluida en ella no 
puede producir la excitación del tejido, mientras que una despolarización en una región en 
cual la región crítica se encuentra estrictamente incluida, puede producir excitación. El 


á 2 . ez Les e 
valor propio ¿ti correspondiente a una región crítica es el valor crítico 


Da 
212 Li 


z 1. La combinación de parámetros geométricos y bio-eléctricos 
4.c 


Mic 


que hacen que una región de influencia hemi-elipsoidal sea una región crítica verifica la 


8-41, 
igualdad: Eu («2 ze , y ESOO 


(61) 


+ . 

Am G 4 do s.m 
Para cada región de influencia despolarizada por acción de un foco ectópico o de un elec- 
trodo externo, que incluya estrictamente a una región crítica, la excitación eléctrica del 
tejido es posible y hay un estado umbral para el sincicio con anisotropía desigual. Hay 
muchas regiones de influencia en esas condiciones, y a cada una le corresponde un patrón 
espacial umbral de despolarización de las membranas excitables. Por tanto, al igual que 
para el sincicio isótropo estudiado en [20], no hay un único estado umbral, sino una 
familia de estados umbral. 


(6.2) Curvas intensidad duración para el modelo con anisotropía desigual. 
Si una órbita atraviesa la barrera umbral (57), proveniente del conjunto de puntos accesi- 


bles situados por encima de dicha barrera, se desencadena un potencial de acción: la esti- 
mulación es exitosa. Si retorna al origen sin alcanzar la barrera umbral no se produce el 
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potencial de acción: la estimulación fracasa. En principio, de acuerdo con este modelo, 
una órbita puede alcanzar la barrera umbral comenzando en un estado accesible tanto por 
aplicación de un pulso catódico (excitación catódica de cierre) como por aplicación de un 
pulso anódico (excitación anódica de apertura). La excitación catódica de apertura, la 
excitación anódica de cierre y la excitación anódica de apertura debido al efecto del mo- 
vimiento de los cátodos virtuales una vez interrumpido el pulso anódico, requieren añadir 
al menos un modo adicional al análisis modal, en forma análoga a como se hizo en las 
referencias [18] y [19] para la estimulación de fibras nerviosas con el propósito de estu- 
diar situaciones análogas aunque bastante más simples que las que se deben considerar en 
el caso de la electro-cardiología. 

El propósito de esta serie de trabajos sobre la estimulación eléctrica del miocardio, es 
estudiar, desde un punto de vista teórico-analítico cualitativo, las dinámicas justo hasta el 
umbral, pero no más allá. Como consecuencia se puede hacer una aproximación lineal a 
la dinámica en el conjunto de decaimiento, en torno al origen, reteniendo la barrera 
umbral, como se hizo y fundamentó en las referencias [14], [17], [18], [19] y [20]: 


dA IB 
Tu, A+ Ki cie) (62 a) ty AB (62 b) 


Introduciendo un vector columna de estados x(+) de componentes A; (1) y B; 0) , Junto con 
la distancia p entre el estado de reposo (el origen del plano de estados) y la barrera um- 
bral, la condición necesaria y suficiente para que en algún instante 1, el estado x(t,) pase 
por la barrera umbral es que 1 x(t,)= p , donde n es un vector unitario columna en el 


plano de estados que es ortogonal a la barrera umbral y apunta hacia el interior del con- 


junto amplificador, y 1? es el vector unitario traspuesto (versor fila). Ver Figura 6. 


B, 


JS Barrera umbral 


Figura 6. Esquema de la construcción geométrica que se utiliza para deducir la condición que 
debe cumplir un pulso estimulador justo umbral. 


La integración del sistema de ecuaciones diferenciales (62 a) y (62 b) junto con la aplica- 
ción de la condición necesaria y suficiente para que un pulso estimulador que comienza 
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en el instante 1 =0sea justo umbral max hi"z(0))= p , siguiendo el procedimiento desa- 
tel0,+00) 

rrollado en la referencia [19], después de varios pasos de cálculo que son formalmente 

iguales a los detallados en la mencionada referencia, permite deducir la siguiente condi- 

ción para un pulso justo umbral aplicado a través de un electrodo de marcapaso: 


max Ú Cl) de) da =01.0 (63) 


te|0,+20) | O 
La condición (63) es un ejemplo de la funcional de excitación introducida por el autor en 
el marco de la cardio-estimulación eléctrica crónica [12] [13] y extendida posteriormente 
a la estimulación de fibras nerviosas y musculares esqueléticas [14] [18] [19]. 
Aquí la función de respuesta impulsiva Gl) viene dada por: 


E)- gt E gota (64) 
K] =K2 K] —K2 
Por definición: k, = - diem?) (65a) k, = a (65 b) 
M ae 


s A 2 
En general, como 7,, es un orden de magnitud mayor que 7,, , mientras que 1+ 4; es 
mayor que y , cabe esperar que x, sea un orden de magnitud mayor que k,. 


La figura 7, tomada de la referencia [19] (pero aplicable también en el caso que estamos 
considerando desde un punto de vista cualitativo), presenta las propiedades más importan- 
tes de la función de respuesta impulsiva. 


1.0 


Gt) 


sd AAA 


u m tiempo 


Figura 7. Gráfica de la función de respuesta impulsiva donde aparecen el tiempo útil hasta la 
reobase í,, y el instante /,, donde la función presenta su mínimo. 
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Para un pulso rectangular de corriente catódica aplicado durante el intervalo [o, t sl: de la 


ecuación (63) se desprende de una fórmula para la correspondiente curva intensidad- 
duración: 


li. (tp)= Beti Up = ta ) (66 a) lil..x = li]. (e, E %o (tp > l, ) (66 b) 


] Cerda ptas 


De la ecuación (66 a) y del signo positivo de Gl) mientras fp <f,, se deduce que la 
amplitud umbral li]. (1,) decrece monótonamente a medida que la duración del pulso se 
incrementa. De la ecuación (66 b) se deduce que la amplitud umbral permanece constante 
en su valor mínimo [:l.. p (reobase catódica) una vez que la duración del pulso supera el 


valor 1, (de ahí la designación de tiempo útil hasta la reobase). El tiempo útil hasta la 


reobase de la Figura 7 se puede estimar a partir de la fórmula %.: Pra IAE y el 


instante en el que se produce el mínimo de Gl) se puede estimar mediante ; =2,, 
El parámetro Oy, y se puede interpretar como el límite de la carga que debe entregar un 


pulso catódico justo umbral de duración tendiente a cero: Op o = limltp «Li ls (17) 
] tp YO 


En el caso del modelo de miocardio con anisotropía desigual, estimulado mediante un 


A 


electrodo de marcapaso, la carga umbral límite viene dada por:  Q,.,= Tm" (67) 


1 


La constante de tiempo para la estimulación catódica f, . = Re y el denominado pará- 


am 


también se expresan mediante fórmulas analíticas en 
i 
c,R 


función de los parámetros del modelo matemático del umbral. En las referencias [14] y 
[19] se presenta un análisis matemático de las curvas amplitud umbral versus duración del 
pulso anódico, para la excitación anódica de apertura debida al efecto de la variable de 
recuperación, que es formalmente aplicable también en este caso. Pero como en el caso 
del miocardio cabe esperar que x, sea al menos un orden de magnitud mayor que x,, el 


metro sigma catódico o, = 


umbral de estimulación anódica de apertura por el mecanismo que involucra solamente la 
variable de recuperación de la membrana origina amplitudes umbral que superan amplia- 
mente las originadas por el mecanismo de excitación anódica de apertura por desplaza- 
miento de los cátodos virtuales. Este último mecanismo depende del grado elevado de 
anisotropía desigual que presenta el miocardio ventricular considerado en el presente 


trabajo. 


=k; +1 


Si se desprecia x, respecto de k;, Glt)= e En este caso las fórmulas (66 a) y (66 b) 


se reduce a una sola: li]. (c,)= e (68) 


Le 


l-e 
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La constante de tiempo viene dada por la fórmula: li ¿7 - M z (69) 
+A] 
2 
á a k (+ 22) Lu 
La reobase correspondiente verifica: li]. po E (70) 


1 
En principio las fórmulas (68), (69) y (70) se pueden aplicar a un miocardio heterogéneo 
con anisotropía desigual, si se calculan los parámetros que allí aparecen. Para un modelo 
homogéneo y anisótropo con propiedades bio-eléctricas iguales en las direcciones per- 
pendiculares a las fibras, teniendo en cuenta la expresión (59) para u de , las fórmulas (69) 
y (70) se reducen a: 


a yd (71) 
e E a. Cs Ja 
8 Xm j G; ; 
E 7” Ry (a Cua c.) ' e A 
[po A (2) 
d?K, 


(7) Conclusiones 


- El concepto de región crítica fue extendido a un sincicio heterogéneo y anisótropo y se 
obtuvo la fórmula (61) para estimar sus dimensiones en un caso de sincicio homogéneo 
con anisotropía desigual del tipo que se observa en el miocardio ventricular. Puede apli- 
carse para interpretar la región liminal hallada experimentalmente por Lindemans [7] y la 
fórmula (61) puede utilizarse para predecir sus dimensiones. Asimismo, se confirma la 
predicción de una mayor vulnerabilidad del miocardio respecto de posibles focos ectópi- 
cos que despolaricen las membranas en una región adyacente a la frontera, que frente a 
focos ectópicos de la misma intensidad que despolaricen las membranas en una región de 
igual volumen alejada de la frontera. 


-Sustitúyase d por y: (r, + e) en la fórmula (71) para la constante de tiempo de la curva 


intensidad-duración por excitación catódica de cierre, donde r, es el denominado radio 


eléctrico del electrodo [13] [14] [15] [16] , e es una medida del espesor del casquete de 
tejidos no excitables interpuestos entre la superficie del electrodo de marcapaso y el mio- 
cardio excitable y y es un número próximo a uno que es mayor en la estimulación bipolar 


respecto de la unipolar. El resultado que se obtiene explica la mayor parte de los resulta- 
dos experimentales sobre las constantes de tiempo para la estimulación catódica de cierre 
[10] [13] [14] [15] y predice algunos resultados nuevos que tendrían que ser contrastados 
contra los datos obtenidos en experimentos específicos. 
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Ar disminuye rápidamente (res- 
d?.K, 
pecto de la escala en la que varía d?) hacia un límite independiente de d . Como conse- 
cuencia, la reobase, en función de d , se expresa como la suma de un término práctica- 


-Téngase en cuenta que, a medida que d aumenta, 


mente independiente de d y un término prácticamente proporcional a d?. Sustitúyase en 
la fórmula (72) d por xr +e). Con esta sustitución la fórmula (72) explica los resul- 


tados experimentales sobre la dependencia de la reobase, para la excitación catódica de 
cierre, del radio del electrodo, el espesor del casquete de tejidos no excitables y la moda- 
lidad (uni o bipolar) de estimulación cardíaca [13] [14] [15]. 

-Se confirma lo expresado en la referencia [20] sobre el concepto de región de influencia 
de un electrodo externo sobre el miocardio, sobre la familia de estados umbral para un 
mismo corazón, en un mismo estado fisiológico o fisiopatológico (la familia de estados 
umbral depende de la forma, dimensiones y distancia del electrodo al miocardio) y la 
distinción entre un concepto global del umbral de excitabilidad del miocardio y un con- 
cepto local de ese umbral. No obstante, la región de influencia cuando se tiene en cuenta 
la anisotropía desigual depende fundamentalmente del factor de forma de volumen, que 
no existe en un modelo homogéneo e isótropo como el estudiado en la referencia [20]. 
-Resta por analizar la dinámica del umbral añadiendo amplitudes de modo de orden supe- 
rior, para poder investigar las excitaciones anódica de cierre, catódica de apertura y 
anódica de apertura sobre el miocardio ventricular. Todo esto será motivo para otro traba- 
jo. 
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Resumen: En esta tercera parte de la serie se utilizan algunos de los resultados obtenidos en los dos 
artículos precedentes para analizar el sistema electrodo-miocardio y discutir, en el marco de este 
sistema, la influencia del tamaño, la forma, el material y la textura de los electrodos de marcapaso 
sobre los umbrales de captura tal como se determinan en la práctica clínica y se utilizan en la fija- 
ción de márgenes de seguridad aplicables en la estimulación crónica. Se indica brevemente cómo 
los resultados de este análisis pueden utilizarse para el diseño de electrodos de marcapaso y se 
presenta un ejemplo: el electrodo clotoide. 


Palabras clave: Marcapasos cardíacos. Sistema electrodo-miocardio. Diseño de electrodos de 
marcapaso. Umbrales de excitabilidad globales y locales. Familias de estados umbral. Regiones 
liminales. Estimulación catódica de cierre. Curvas intensidad-duración. Electrodos cóncavos. Elec- 
trodo clotoide. 


Abstract: In this third part of the series some of the results obtained in the previous papers are used 
to analyze the electrode-myocardium system and to discuss, in the framework of this system, the 
influence of the size, shape, material and texture of pacemaker electrodes on capture thresholds 
such as there are determined in clinical practice and used to establish safety margins in chronic 
cardiac pacing. Some brief indications about how the results of this analysis can be applied in 
pacemaker electrode design and one example is given: the clothoidal electrode. 


Key words: Cardiac pacemakers. Electrode-myocardium system. Pacemaker electrodes design. 
Global and local excitability thresholds. Families of threshold states. Liminal regions. Cathodic 
make stimulation. Strength-duration curves. Concave electrodes. Clothoidal electrode. 


(1) Introducción 


La introducción del marcapaso cardíaco artificial para tratar trastornos del ritmo cardíaco 
fue una de las innovaciones médicas más importantes del siglo pasado. 

Cuando el marcapaso natural del corazón no funciona bien y se produce una arritmia, los 
impulsos eléctricos del marcapaso artificial pueden restituir el ritmo cardíaco a valores 
adecuados. Como la estimulación eléctrica del corazón presenta un umbral, si se fija una 
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duración para el impulso, su amplitud debe superar un valor mínimo para que el equipo 
pueda capturar el miocardio, provocando un potencial de acción que se propague a toda la 
masa de músculo y pueda desencadenar una contracción ventricular adecuada para man- 
tener la función de la bomba cardíaca. 

Los marcapasos constan de un generador de impulsos eléctricos, un cable (o cables) en el 
interior de un catéter de un material aislante y un electrodo (o electrodos) en el extremo 
del cable (o cables). 

Cuando el electrodo en el extremo es único, el circuito eléctrico se cierra, a través del 
conductor de volumen de los tejidos, entre ese electrodo que actúa como cátodo y un 
ánodo formado por parte de la superficie externa del generador de pulsos. La correspon- 
diente modalidad de estimulación se denomina unipolar. Cuando el catéter presenta, ade- 
más del electrodo en su extremo (cátodo), otro electrodo próximo al primero (ánodo), el 
circuito eléctrico a través de los tejidos se cierra entre esos dos electrodos. En este caso, la 
superficie externa del generador de pulsos se encuentra completamente aislada desde el 
punto de vista eléctrico. La modalidad de estimulación se denomina bipolar [1], [2]. 


Los primeros marcapasos artificiales relativamente exitosos, empleados en pacientes con 
la intención de controlar el ritmo cardíaco en forma crónica, fueron construidos durante la 
primera mitad del siglo pasado. Tenían la unidad generadora de los trenes de impulsos 
eléctricos ubicada fuera del organismo. Esto limitaba mucho su practicidad y aumentaba 
el riesgo de infección, puesto que era necesario atravesar la piel con un cable para alcan- 
zar el corazón con un electrodo. Los intentos por implementar un marcapasos totalmente 
incorporable en el interior de un organismo culminan en 1958 con la primera implanta- 
ción clínica de un marcapasos artificial interno, diseñado por Rune Elmgqyvist, en Suecia 
[2]. A comienzos de 1960, Orestes Fiandra y Roberto Rubio implantaron a un paciente, 
en Uruguay, una versión mejorada del marcapasos sueco, también debida a Elmqyvist. Se 
considera el primer implante de un marcapasos en América [3]. 


Durante los primeros 40 años del desarrollo de los marcapasos, tanto los ingenieros como 
los médicos tuvieron bastante éxito en definir y concretar las necesarias mejoras en el 
hardware y en el software, muchas veces por ensayo y error, sin una comprensión a fondo 
de los mecanismos subyacentes a la estimulación del músculo cardíaco [4]. 

Para avanzar en el diseño, construcción y empleo de los marcapasos artificiales se requie- 
re una comprensión más profunda de los mecanismos subyacentes a la emergencia y pro- 
pagación de los potenciales de acción en el miocardio, en las condiciones específicas 
debidas a la interacción entre los electrodos y el tejido miocárdico. La mayor parte de los 
problemas de investigación todavía abiertos involucran directa o indirectamente a los 
electrodos de marcapaso: la interfase entre la máquina y el corazón [1], [4]-[7]. 


Como se dijo en los dos artículos precedentes [4], [7], el propósito principal de esta serie 
es contribuir a una mejor comprensión de los problemas del umbral de estimulación eléc- 
trica del miocardio a partir del análisis de un modelo matemático de la interacción entre el 
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electrodo de marcapaso y el tejido cardíaco excitable. 

En esta tercera parte se utilizan algunos de los resultados obtenidos en los dos artículos 
mencionados para analizar el sistema electrodo-miocardio y discutir, en el marco de este 
sistema, la influencia del tamaño, la forma, el material y la textura de los electrodos de 
marcapaso sobre los umbrales de captura tal como se determinan en la práctica clínica y 
se utilizan en la fijación de márgenes de seguridad aplicables en la estimulación crónica. 
Se indica brevemente cómo los resultados de este análisis pueden utilizarse para el diseño 
de electrodos de marcapaso y se presenta un ejemplo: el electrodo clotoide. 


Para comenzar conviene revisar algunos aspectos básicos. La Figura 1 muestra un esque- 
ma de un marcapaso con un único electrodo próximo al miocardio (electrodo unipolar), 
ubicado en el interior del ventrículo derecho, apoyado en el endocardio (estimulación 
endocárdica). 


Catéter 


Generador de 
pulsos 


yr "Y 
42 
". 


NES 
y Cabeza del 


electrodo 


Figura 1 Esquema de un marcapaso con electrodo unipolar. La cabeza del electrodo actúa 
como cátodo y el generador de pulsos implantado actúa como ánodo. 


La Figura 2 presenta un esquema del catéter con el cable conductor en su interior, la 
cabeza del electrodo junto con elementos de fijación pasiva al miocardio en un extremo y 
el conector con el generador de pulsos en el otro extremo. 
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Cabeza Cable Conductor  Aislante Conector 


Figura 2 Esquema de un catéter con la cabeza de electrodo en uno de sus extremos y un co- 
nector para el generador de pulsos en el otro. Modificado de [2]. 


El catéter con el cable se puede hacer ingresar a través de la vena subclavia izquierda, 
entrando en la aurícula derecha a través de la vena cava superior, atravesando la corres- 
pondiente válvula aurículo-ventricular y alojando su extremo en el vértice anfractuoso del 
ventrículo derecho, donde los elementos de fijación del extremo del catéter contribuyen a 
mantenerlo en posición. 

El par de electrodos (el electrodo unipolar en el extremo del catéter y el electrodo situado 
en la superficie externa del generador de pulsos) se emplea tanto para estimular el cora- 
zón como para medir la actividad eléctrica espontánea del miocardio. 

En algunos casos se estimulan dos o más cámaras cardíacas utilizando varios catéteres en 
simultáneo [2]. 


En la estimulación crónica el generador de pulsos produce un tren de pulsos de voltaje 
controlado entre el ánodo y el cátodo, con una separación temporal entre pulsos y una 
duración de pulso adecuadas para el caso. 

Si se fija una duración para el pulso estimulador, su amplitud debe superar un valor mí- 
nimo (umbral) para que el equipo pueda capturar el miocardio, provocando un potencial 
de acción que se propague a toda la masa de músculo. El umbral se suele determinar apli- 
cando el pulso fuera de los períodos refractarios absoluto y relativo (umbral diastólico 
tardío). A menor duración del pulso, mayor resulta la amplitud umbral. 


Si U,, representa la amplitud umbral de un pulso de voltaje controlado y £ pes su dura- 


ción, la dependencia entre ambos se describe a menudo en la práctica clínica mediante la 
fórmula hiperbólica de Weiss-Lapicque: 


O E — (1) 
P 
En esta expresión el umbral se caracteriza mediante dos parámetros: U y, que es la am- 
plitud reobásica (la amplitud umbral mínima, por debajo de la cual por larga que sea la 
duración del pulso estimulador, la captura ya no es posible) y í,, y que es la cronaxia 
(duración del pulso de voltaje para la cual la amplitud umbral es el doble de la reobase). 


De acuerdo con la fórmula (1), la amplitud umbral disminuye y se aproxima a la reobási- 
ca a medida que la duración del pulso aumenta tendiendo a infinito. En la práctica la reo- 


Rev. Soc. Cient. Paraguay. Tercera Época. 2016; 21(2) ISSN 0379-9123 


base de voltaje se alcanza para una duración finita (duración útil hasta la reobase de volta- 
je) que a menudo está próxima a los 2 ms. Las cronaxias toman valores de décimas de ms, 
generalmente valores comprendidos entre 0.2 y 0.5 ms. Las reobases de voltaje varían 
más, pero en general toman valores de voltios. 


También se determinan los umbrales de pulsos de corriente controlada para diferentes 
duraciones de pulso. En la práctica clínica los resultados de las mediciones del umbral de 
corriente se suelen expresar también recurriendo a una fórmula de Weiss-Lapicque: 


t 
L,=Lg | 1+ e Q) 


Pp 


En esta fórmula / y, es la reobase de corriente (que en la práctica también se alcanza para 


una duración finita, la duración útil hasta la reobase de corriente) y £,, ,es la correspon- 


diente cronaxia de corriente. Para un mismo sistema electrodo-miocardio, cuando se las 
determina en el mismo momento, por lo general la cronaxia de voltaje resulta algo menor 
(o a lo sumo igual) que la cronaxia de corriente [8]. 

En trabajos que describen los resultados de las determinaciones de curvas amplitud um- 
bral versus duración para pulsos rectangulares de corriente controlada, más orientados 
hacia la electrofisiología cardíaca, se utiliza a veces la fórmula de Lapicque-Hill. El um- 
bral de corriente se caracteriza mediante la misma reobase y la denominada constante de 
a través de la expresión: Jia A 6) 


u t 


tiempo £, ,, 


La constante de tiempo es proporcional a la cronaxia: to, = (In 2)- bei (4) 


Para una misma duración del pulso, los resultados experimentales aportan evidencia acer- 
ca de que los umbrales dependen, entre otras cosas de los siguientes factores [2], [6], [9]- 
[12]: 


(1) La ubicación del electrodo en el interior de la cámara cardíaca estimulada. 


(Q) El estado del miocardio excitable, influido por metabolitos, señales endócrinas y neu- 
rales, globales y locales. Los ritmos circadianos, el sueño, el consumo de alimentos, la 
actividad física, la fiebre, las drogas, las emociones y otros factores inciden sobre los 
valores del umbral. 


(3) La reacción inflamatoria desencadenada por la presencia del electrodo, que termina 
con la formación de un casquete de tejido fibroso no excitable en torno a la cabeza del 


electrodo. 


(4) La mayor o menor rigidez del catéter. 
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(5) Las propiedades del electrodo que se encuentra en su extremo (tamaño, forma, mate- 
rial y textura). 


(6) La distancia entre el electrodo y el tejido excitable. 


(7) El procedimiento utilizado para determinar el umbral. 

Cuando el umbral se mide comenzando por pulsos supra-umbrales cuya amplitud se dis- 
minuye hasta que el marcapaso pierde la captura del miocardio, los umbrales así determi- 
nados resultan inferiores respecto de los valores que se obtienen cuando el umbral se 
mide comenzando por pulsos sub-umbrales cuya amplitud se va incrementando hasta que 
el equipo logra capturar la actividad eléctrica del miocardio [13]. 


En consecuencia, las curvas amplitud umbral versus duración para pulsos de corriente y 
de voltaje controlado son una propiedad del sistema electrodo-miocardio considerado 
como un todo. Describen lo que puede denominarse excitabilidad global de ese sistema. 
La excitabilidad global debe distinguirse de la excitabilidad local. Esta última es la exci- 
tabilidad definida a nivel de las membranas de las fibras miocárdicas donde se produce el 
patrón de polarización originado por el pulso estimulador. 


La excitabilidad global se relaciona con la excitabilidad local a través de ciertas regiones 
críticas de células miocárdicas, cuyas membranas deben despolarizarse por encima del 
umbral uniforme de membrana para que se pueda originar un potencial de acción capaz 
de propagarse al resto de la masa del miocardio [4], [7], [14], [15]. 

Para discutir (teniendo en cuenta la distinción entre excitabilidad global y excitabilidad 
local) la influencia del tamaño, la forma, el material y la textura de los electrodos de mar- 
capaso sobre los umbrales de captura, es conveniente comenzar por una descripción fun- 
cional del sistema electrodo-miocardio. 


(2) Descripción funcional del sistema electrodo-miocardio. 


Los electrodos y los tejidos biológicos del paciente forman un sistema cuyas propiedades 
eléctricas, electroquímicas, mecánicas y biológicas se deben tener en cuenta en el proceso 
de diseño de los electrodos de marcapaso. 

Si se enfatiza el flujo de señales a través del sistema, se pueden distinguir cuatro subsis- 
temas parcialmente superpuestos en el espacio. Ver Figura 3. 
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AREA BIOELECTRICAPASIVA 
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AREA BIOELECTRICA ACTIVA 


Figura 3: Esquema para una descripción funcional del sistema electrodo-miocardio [28]. 


(A)Un subsistema eléctrico pasivo, que denominamos área bioeléctrica pasiva. Está 
compuesta por cuatro módulos conectados en serie: 

(1) El catéter con el cable y el electrodo activo (cátodo) en su extremo proximal al mio- 
cardio (a lo que se añade el cable del ánodo y el ánodo si la estimulación es bipolar). 

(Q) La interfase entre la superficie del electrodo activo y los tejidos biológicos. 

(3) Un casquete de tejidos no excitables que se interpone entre la superficie del cátodo y 
el miocardio excitable, en la fase crónica de la estimulación consiste en una cápsula de 
tejido fibroso. 

(4) El conductor de volumen formado por el resto de los tejidos, hasta el electrodo indife- 
rente. 

En el área bioeléctrica pasiva se pueden distinguir dos señales de entrada y una señal de 
salida. Una de las señales de entrada es eléctrica: es un pulso de corriente que atraviesa 
cada uno de los cuatro módulos situados entre el electrodo activo y el electrodo indiferen- 
te. La otra señal de entrada es cinemática, proviene del cuarto subsistema y será conside- 
rada luego de describir dicho subsistema. La señal de salida es una distribución espacial 
de corriente eléctrica que conecta este subsistema con el subsistema siguiente. 


(B)Un subsistema eléctrico activo, que denominamos área bioeléctrica activa. Está 
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compuesto por tres módulos, con los dos primeros conectados en serie: 
(1) El tejido cardíaco excitable situado próximo al cátodo, blanco de la estimulación. 
(2) El resto del miocardio. 
(3) La jerarquía de marcapasos naturales. 
La señal eléctrica de salida desde el área bioeléctrica pasiva (el patrón de distribución de 
la densidad de corriente eléctrica en el conductor de volumen de los tejidos) es la señal de 
entrada al área bioeléctrica activa. En las membranas de las células miocárdicas próximas 
al electrodo la distribución de corrientes produce un patrón de polarización característico 
de la anisotropía desigual del tejido, considerada en uno de los artículos previos [39]. En 
el área bioeléctrica activa se pueden distinguir dos señales de salida: un patrón de ECG 
cuyas características permiten constatar si la captura del miocardio se ha producido o no, 
y una señal físico-química (potencial de acción) que conecta este subsistema con el si- 
guiente. 


(C)Un subsistema mecánico activo o área biomecánica activa. Formado por el músculo 
cardíaco contráctil, recibe la señal fisico-química del área bioeléctrica activa y se con- 
trae produciendo una señal mecánica consistente en un campo variable de esfuerzos 
que lo conecta con el subsistema siguiente. 


(D)Un subsistema mecánico pasivo o área biomecánica pasiva. Incluye los tejidos que 
rodean la cabeza del electrodo activo y el catéter, considerados desde el punto de vista 
mecánico. 

Recibe la señal mecánica proveniente del área biomecánica activa y produce una señal de 

salida cinemática. Esta señal cinemática entra en el área bioeléctrica pasiva donde se 

evidencia en una modificación del espesor de la cápsula de tejido no excitable interpuesto 
entre la superficie del electrodo activo y el miocardio excitable blanco de la estimulación. 

Esto hace variar la distancia electrodo-miocardio que aparece sumada al radio eléctrico 

del electrodo en las fórmulas para las reobases y las constantes de tiempo obtenidas en los 

dos artículos previos [4], [7]. 


Se cierra así un circuito de transferencia de señales, que comenzando con la entrada del 
pulso de corriente en el área bioeléctrica pasiva termina con una señal cinemática de en- 
trada a esta misma área, proveniente del área biomecánica pasiva. 

Este bucle de retroalimentación complica la construcción de un modelo matemático del 
sistema electrodo-miocardio durante la estimulación por electrodos de marcapaso, debido 
a que la distancia entre el electrodo y el miocardio excitable puede variar de un pulso a 
otro, dependiendo del grado de complacencia del catéter, de la frecuencia del tren de 
pulsos estimuladores, de la amplitud de la contracción del músculo cardíaco, de la ubica- 
ción y anclaje del electrodo, entre otros parámetros. 


No obstante es posible abrir el lazo de retroalimentación considerando la distancia elec- 
trodo-miocardio excitable como una variable aleatoria, que toma un valor constante du- 


Rev. Soc. Cient. Paraguay. Tercera Época. 2016; 21(2) ISSN 0379-9123 


rante la duración de un mismo pulso estimulador, aunque ese valor pueda modificarse de 
un pulso a otro. Esta constancia para un mismo pulso se debe a que los pulsos poseen una 
duración del orden de las milésimas de segundo, mientras que las fases del ciclo cardíaco 
y los fenómenos mecánicos en el sistema electrodo-miocardio se desarrollan en escalas 
del orden de las décimas de segundo. 

Una consideración detallada de las escalas de tiempo que se pueden distinguir en la esti- 
mulación cardíaca crónica y que conviene tener en cuenta durante la construcción de 
modelos matemáticos de este proceso (microsegundos, milisegundos, décimas de segundo 
y horas) puede hallarse en las referencias [16], [17]. 


(3) La funcional de excitación y los pulsos umbral. 


A partir de la formulación de las ecuaciones acopladas que describen el flujo de corriente 
eléctrica en el conductor de volumen de los tejidos y la polarización de las membranas 
excitables basada en el modelo de bidominio [18], de la región de influencia del electrodo 
sobre el miocardio excitable [15], de los resultados del análisis modal no lineal del umbral 
aplicado a las ecuaciones de bidominio (utilizando la región de influencia) y de la identi- 
ficación de un modo dominante [4,7], se obtiene una condición que deben cumplir los 
pulsos de corriente eléctrica ie ) introducidos en el área bioeléctrica pasiva para poder 
capturar la actividad eléctrica del miocardio justo en el umbral. Esa condición se formula 
en términos de una correspondencia que a cada función ie ) le hace corresponder un 
escalar (una carga umbral límite): una funcional a la que el autor propuso denominar 
funcional de excitación [4], [7], [14], [15], [19], 20]. 

Si esta expresión se aplica a un tren de pulsos lo bastante separados como para que la 
influencia del pulso precedente no sea significativa, es suficiente considerar un único 
pulso aislado para analizar los fenómenos relacionados con el umbral. Por lo general, éste 
es el caso en la estimulación mediante electrodos de marcapaso. El pulso de corriente 
puede ser tanto un pulso de corriente controlada como el pulso de corriente asociado con 
un pulso de voltaje controlado. 


El criterio de umbral basado en la funcional de excitación vincula la entrada al área bio- 
eléctrica activa con una de las dos salidas de dicha área, la del tipo todo o nada, según que 
la captura del miocardio haya tenido o no éxito. Como se dijo en la sección 2 precedente, 
la respuesta de tipo todo o nada se pone en evidencia a través del patrón de ECG [5]. 

La funcional de excitación depende de los parámetros del sistema electrodo-miocardio 
como se explicó en la referencia [7]. 


En general, tanto en relación con la disminución del umbral como en relación con la 
energía disipada en los tejidos, en igualdad de las demás condiciones las formas de pulso 
de corriente que presentan una fase creciente hasta el final del pulso (por ejemplo, una 
forma trapezoidal creciente) son mejores que las formas de pulso de corriente que presen- 
tan fases decrecientes (por ejemplo, una forma trapezoidal decreciente) [11], [15]. 
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(4) Análisis de los umbrales de captura mediante pulsos de corriente controlada. 


Para un pulso rectangular de corriente catódica, después de varias simplificaciones que se 
desarrollan en la referencia [7], a partir del criterio de umbral basado en la funcional de 


excitación se desprende esta fórmula para la curva intensidad-duración, donde / , repre- 


senta la amplitud umbral de corriente y f, representa la duración del pulso: 


1. (ty) 2 (5) 
l-e * 
2 
h (e 
La constante de tiempo viene dada por [7]: 2 úl d, (6) 
SET ATAR 
1 de 
d, 
La reobase viene dada por [7]: Lp (a; A a As 3 (7) 
K,:d 


En estas fórmulas 7,, es una constante de tiempo característica de la polarización uni- 
forme de la membrana de la fibra miocárdica, d = Y * (r, + e) es una distancia formada 


a partir del radio eléctrico », del electrodo [15, 21], de una medida del espesor e del 


casquete de tejidos no excitables, y de un factor sin dimensiones Y que toma valores 
próximos a uno!. Por otra parte, d o €s la distancia para la cual la constante de tiempo 


£. , es la mitad de la constante de tiempo de membrana. 


De los resultados obtenidos en la referencia [7] para un modelo de miocardio homogéneo 
con anisotropía desigual, se desprende que d y es función de los parámetros que caracte- 


rizan las fracciones de volumen y las conductividades de los dominios intracelular y ex- 


| Se puede estimar y 4 para la estimulación monopolar con cátodo convexo [19], [20]. Es mayor 
3 

en la estimulación bipolar respecto de la unipolar, y en general cabe esperar que aumente cuando las 

líneas de corriente en los tejidos tiendan a separarse respecto del eje del cátodo [14], [15]. 
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tracelular, de la fracción de área YX, de membranas excitables por unidad de volumen 


del bidominio cardíaco y de la resistencia R,, a la polarización de la unidad de membra- 
na excitable en las proximidades del potencial de membrana de reposo?, Una estimación 


rápida del orden de magnitud de d, se puede obtener asumiendo isotropía. En ese caso 


2 
45-71" Ry -C= E «7-4, Siendo A, la constante de espacio del cable tridi- 
Xu 


mensional equivalente al bidominio cardíaco. Poniendo 2,, =0,6mm ([14], [15]) se ob- 
tiene d, =2,23mm. 


Del análisis modal desarrollado en la referencia [7] se desprende que 4,, es un valor 


Lu 
umbral de un modo dominante y K, describe el efecto de la geometría de la distribución 


de la densidad de corriente sobre la dinámica del modo dominante. 


El producto de la reobase de corriente por la correspondiente constante de tiempo es la 


Tu Ars 
(8) 
K 


1 


carga umbral límite: O IA 


De la fórmula (6) se desprende que la constante de tiempo f; , debería aumentar cuando 


d= Ye (r, + e) aumenta. Entonces, igualdad de las demás condiciones, la constante de 


tiempo debería aumentar con el radio del electrodo y con el espesor del casquete de teji- 
dos no excitables interpuesto entre el electrodo y el miocardio. También cabe esperar que 
sea algo mayor en el caso de estimulación bipolar, puesto que, en igualdad de las demás 
condiciones, y es mayor que en el caso unipolar. 


A medida que d aumenta, Ar disminuye rápidamente (respecto de la escala en la 
d*-K, 

que varía d?) hacia un límite independiente de d [7]. Como consecuencia, la reobase, en 

función de d , se expresa como el producto de un término prácticamente independiente de 

d por un término proporcional a d¿+d”. 


2 Para un modelo de bidominio homogéneo y anisótropo, con anisotropía desigual, se obtiene la 
siguiente expresión para la distancia característica, cuando se supone que las fibras se disponen 
paralelas al plano que representa la frontera del miocardio próxima al electrodo y se asume un 
comportamiento eléctrico isótropo para direcciones perpendiculares al eje de las fibras: 


2 
don Es Lu (e, qa) 
Xa G, 4 


El significado de los parámetros eléctricos C,,C,,G,»G, se explica en las referencias [4], [7]. 
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Entonces, de la fórmula (7) se desprende que, en igualdad de las demás condiciones, la 
reobase debería aumentar con el radio del electrodo, con el espesor del casquete de tejidos 
no excitables interpuesto entre el electrodo y el miocardio y con el cambio en la modali- 
dad de estimulación de monopolar a bipolar. 


En general, para los electrodos convexos convencionales, la evidencia experimental con- 
firma estas predicciones, tanto en el caso de la constante de tiempo como en el caso de la 
reobase [1], [6], [8], [9], [11], [12], [22], [23]. 


Al final de la sección 2 se dijo que la distancia electrodo-miocardio se podía considerar 
como una variable aleatoria: pese a que es posible suponer que esta distancia variable e 


toma un valor constante durante la duración de un mismo pulso estimulador, ese valor 
pueda modificarse aleatoriamente de un pulso a otro. 


Es conveniente descomponer € en la suma de un valor promedio y una fluctuación alea- 
toria en torno a ese valor medio: é=8 + É 

Representando por E [ ] la operación de tomar el valor medio (valor esperado o esperan- 
za matemática) de una variable aleatoria, se tiene E [e]- 0 puesto que E [e] =68, 


Siendo aleatoria la distancia entre la cabeza del electrodo y el miocardio excitable, la 
reobase 1 y, , que depende del cuadrado de esa distancia, también es una variable aleato- 


ria, proporcional a dí +y? ol, +e+ ¿J ; 


Como consecuencia, el valor medio de la reobase es proporcional a 


di + y? Ello +e+ ¿| zz di + y? -(7, +e) +Elé] : en igualdad de las demás 


NN A 22 , 
condiciones, si E E ] aumenta, el valor esperado de la reobase debería aumentar. 


Ahora bien, el valor esperado E lE 2] del cuadrado de la fluctuación en la distancia elec- 


trodo-miocardio depende del comportamiento del área biomecánica pasiva, que es un 
sistema mecánico sobre-amortiguado y forzado por las contracciones del miocardio [28]. 


Al ser un sistema sobre-amortiguado”, la amplitud de la oscilación forzada disminuye 

cuando la frecuencia de la fuerza exterior (en este caso debida a la contracción cardíaca) 
£2 y ao ; 

aumenta. Entonces E E debe decrecer cuando la frecuencia de contracción del mio- 


cardio aumenta. 
Cuando se mide el umbral de estimulación comenzando con pulsos supra-umbrales cuya 


3 Sobre oscilaciones sobre-amortiguadas puede consultarse cualquier buen texto de física general de 
nivel universitario. 
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amplitud se va disminuyendo, se deben emplear trenes de pulsos repetidos a frecuencias 
que son mayores que las frecuencias de contracción espontánea del miocardio a la que se 
determinan los umbrales comenzando con pulsos sub-umbrales cuya amplitud se va au- 


A £2 ; 
mentando. Así pues cabe esperar que E E sea menor cuando el umbral se determina 
comenzando con pulsos supra-umbrales de duración fija y amplitud variable, respecto de 


E lE 2] cuando el umbral se determina comenzando por pulsos sub-umbrales de la misma 
duración y amplitud variable. 


La constante de tiempo f¿ , (o en forma equivalente, la cronaxia) se comporta en forma 


análoga pero comparando la dependencia de la distancia predicha en la fórmula (6) con la 


dependencia predicha en la fórmula (7), cabe esperar que f¿ , varíe mucho menos acusa- 


damente con la distancia en comparación con la reobase 1 y . 


22 
En todo caso, cuando E E ] aumenta, tanto el valor esperado de la reobase como el de 


la constante de tiempo debería aumentar: el umbral de estimulación del miocardio debería 


; 8 O . 
aumentar. Como consecuencia de que E E ] disminuye al aumentar la frecuencia de 


contracción del miocardio, en igualdad de las demás condiciones el umbral de estimula- 
ción debería resultar en promedio menor cuando se mide a partir de pulsos supra- 
umbrales que cuando se mide a partir de pulsos sub-umbrales. 

Además, cuanto mayor es la frecuencia de esos pulsos supra-umbrales, menor debería 
resultar el valor promedio y las fluctuaciones del umbral. Se predice así una componente 
mecánica en el fenómeno de histéresis mencionado en la sección 1. 


Hay una componente puramente electrofisiológica en el efecto de histéresis, relacionada 
con el denominado efecto Wedensky* [13]. Al parecer la componente mecánica analizada 
en el presente trabajo no ha sido tenida en cuenta (con excepción de la referencia [17]). 


Se dispone de alguna evidencia experimental que sugiere que la componente mecánica 
está presente en la histéresis del umbral [11]. Por ejemplo, este fenómeno es más marcado 
para catéteres más rígidos, como los catéteres de Furman que se empleaban cuando el 
marcapaso se utiliza en forma transitoria [11], [17]. 

A medida que la complacencia de los catéteres para implantación definitiva aumenta y los 
medios de fijación de la cabeza del electrodo mejoran, se observa una disminución del 
efecto de histéresis [1], [9]. 


4 En trabajos sobre cardioestimulación eléctrica crónica, muchas veces se utiliza “efecto Wedensky” 
como sinónimo de “efecto de histéresis”. 
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(5) Análisis de los umbrales de captura mediante pulsos de voltaje controlado. 


Los generadores de pulsos de los marcapasos utilizan pulsos de voltaje controlado de 
amplitudes aproximadamente constante (pulsos rectangulares de voltaje). 

Mientras la amplitud del pulso de voltaje entre el cátodo y el ánodo permanece constante, 
el pulso conjugado de corriente eléctrica sufre variaciones que reflejan los cambios en la 
impedancia de carga vista por el generador de pulsos. 

Esta impedancia es independiente de si el miocardio es capturado o no por el marcapaso: 
refleja las propiedades del área bioeléctrica pasiva solamente. Se puede expresar como 
suma de dos términos. 


Uno de ellos es una resistencia R p ue representa el comportamiento óhmico de la com- 
binación en serie de los cables del (o los) electrodo (electrodos) representados por una 
resistencia R, con el conductor de volumen de los tejidos representado por una resisten- 
cia Ry denominada resistencia de los tejidos [24], [17]. 


La resistencia de los tejidos Ry depende del tamaño y la forma del electrodo activo así 

como de la conductividad eléctrica en su vecindad: no es una propiedad de los tejidos en 

sí mismos [17]. 

Si U, es la amplitud de un pulso de voltaje aplicado (no necesariamente umbral) e / (e ) es 

el pulso de corriente conjugado, la resistencia R, = R, +R, se puede determinar así: 
Yo 

t>0 i(e) 


Por otra parte, si 1, es la amplitud de un pulso de corriente introducido a través de los 


R, =lim 


electrodos y si U (e ) es el pulso de voltaje conjugado, R,= lim,_, 10) 
lo 


U,  U() 


En la práctica se suele determinar los cocientes —H o 


(CA 


a los 90 us de iniciado el 
pulso, es decir, para £ =90 us [24]. 
El esquema que aparece en la Figura 4 representa estos pares de señales de entrada y 


salida (pulsos conjugados de corriente y voltaje) en el área bioeléctrica pasiva (compren- 
dida entre el cátodo y el ánodo del marcapaso implantado). 


Rev. Soc. Cient. Paraguay. Tercera Época. 2016; 21(2) ISSN 0379-9123 


Polarization 
potential 


(a) 


(b) 


Figura 4. Formas de pulso conjugadas en el área bioeléctrica pasiva. En (a) se representan dos 


pulsos rectangulares de corriente controlada de amplitud lo junto con los correspondientes 
pulsos conjugados de voltaje U (1) . En (b) se representan dos pulsos rectangulares de voltaje 


controlado de amplitud U, junto con sus pulsos conjugados de corriente (e ). En ambos 


casos los pulsos conjugados presentan una respuesta óhmica seguida de una respuesta diferi- 
da asociada a la polarización de los electrodos de marcapaso. Modificada de [25]. 


Debido a que ie ) disminuye más rápido de lo que U 6 ) aumenta, el valor de R p esti- 
mado para pulsos rectangulares de voltaje controlado resulta un poco mayor que el valor 
de R » estimado para pulsos rectangulares de corriente controlada [24]. Como la resisten- 
cia del cable se puede hallar en forma independiente, una vez conocida R p Se puede de- 


terminar la resistencia de los tejidos. En una primera aproximación esa resistencia se 
puede considerar inversamente proporcional al tamaño del electrodo [10], [17], [21], 26]. 
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El otro término de la impedancia vista por el generador de pulsos del marcapaso, es una 
impedancia que corresponde a la polarización de la capa doble electroquímica que se 
forma en la interfase entre los electrodos y los tejidos biológicos [1], [6], [27]-[29]. 

Desde el punto de vista electroquímico las interfases entre los electrodos de marcapaso y 
los tejidos (o las soluciones fisiológicas) exhiben un comportamiento no lineal. Es necesa- 
rio superar una diferencia de potencial crítica a través de la intefase para que comience a 
fluir una corriente farádica en paralelo con la corriente capacitiva. Una vez superado este 
valor crítico, la relación entre la corriente farádica y la diferencia de potencial que la 
impulsa no es lineal: se parece a lo que se observa conectando en serie dos diodos Zener, 
oponiéndose el uno al otro [17], [28], [29]. Pero para los fines del presente trabajo se 
puede utilizar, en una primera aproximación, un modelo lineal. 

En ese caso las diferencias de potencial U (e) entre el ánodo y el cátodo se pueden escri- 


bir como suma de caídas óhmicas en los cables y los tejidos más las caídas de potencial 
en las capas dobles como se muestra a continuación: 


Ult)=R,-1(0)+ | K((—u)-1(u)-du 03) 


Aquí la función de respuesta al impulso K 6 ) es positiva y caracteriza el comportamiento 


eléctrico de la interfase en una primera aproximación lineal. Si R, representa la resisten- 
00 

cia de la interfase a la corriente directa: J K (s) .dt=R, (9 b) 
0 


Puesto que en las condiciones que se encuentran en electrofisiología cardíaca clínica los 
generadores de pulsos producen pulsos de voltaje tales que por lo general se supera la 
mencionada diferencia de potencial crítica en la interfase, el autor propuso describir el 
comportamiento electroquímico de la interfase mediante una impedancia de Drake [30]. 
Admitiendo esta propuesta, para instantes de tiempo positivos: 


(10 a) 


El parámetro M está comprendido entre 0 y 1. Este parámetro se relaciona con la dimen- 


sión fractal d s de la interfase a través de la acotación: m<3-d; (10 b) 


El parámetro y =| —“% A su vez, el parámetro Ko aparece en la fórmula de la 
d 
impedancia de Drake para la interfase: Lo (o) = —— (10 c) 
(a+ -0) 


Como es usual, (0 representa la frecuencia angular de la corriente alterna y ¡=V—1. 
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T(m)= [ute «e * «du es la función gama de Euler. 


Como el generador de impulsos impone U 6 ), conviene invertir la ecuación (9 a). Se 


A 1 ( 
obtiene: i(s)= z (vw fz( ¿JU a) (11 a) 
Pp 0 
Mi 2 e A ñ R, 
La función Lle) también es positiva y verifica: fz(0)- dt= “GT (11 b) 
o R, +R, 
Cuando K 6 ) se describe mediante la fórmula (10 a) correspondiente a una impedancia 
Tm aa 
de Drake, resulta: L(t) z K(t)- it. : e hr(2-m) Ma) (11 c) 
m-(2-m) 
Cuando se aplica un pulso de voltaje de amplitud constante U,,, la ecuación (11 a) se 
E Ú, 
reduce a la siguiente: i(e) ==. (1 + pls) (12) 
R, +R 
e Ra 


Por definición: plt)= - 14 - | L(u)-du (13) 
p) 0 


De (13) y (11 b) se desprende que la función plt ) es positiva y monótona decreciente, 


R 
desde pl0) = E hasta cero cuando f£ —>00. Entonces de (12) resulta que el pulso de 
P 


corriente 1 (e ) disminuye desde el comienzo hasta el final del pulso de voltaje, tal como 
sugiere la figura 4 (b). De (11 c) se desprende que Ll ) tiende a cero más rápido que 
K 6 ) a medida que pasa el tiempo. Como consecuencia, en la Figura 4 (b) ie ) disminu- 
ye más rápido que lo que aumenta U (e ) en la Figura 4 (a). 

Un pulso rectangular de voltaje será un pulso justo umbral 1, 6 ) si el pulso de corriente 
conjugado de duración f, cumple la condición de umbral establecida empleando el crite- 


rio de la funcional de excitación para un pulso catódico [4, 15]: 


[ale, -1)-i,(0)-di=0,, (14) 
0 


La función Gl: ) depende de algunos de los parámetros del sistema electrodo tejidos 


como se expuso en segundo trabajo de esta serie [39], donde además se puede hallar una 
representación gráfica. 
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Por otra parte, si U, es la amplitud umbral de un pulso de voltaje de duración £ p> Je la 


fórmula (12) se desprende que: E 6 ) = E : (1 + pl ) (15) 
d 


De (14) y (15) resulta la siguiente fórmula para la relación entre la amplitud de voltaje y 
la duración de los pulsos catódicos: 


0.) Eo tR) Qs (16) 


[6(0).ar+ fat, —1t)-p(t)-dt 


Si se aplican las mismas simplificaciones que se mencionaron en la sección 4, para valo- 


t, 


res positivos del tiempo resulta: Gl: ) =8 (17) 
O t, e —1) 
Teniendo en cuenta que /,, =% y que A fe ts -plt)- dt=0, sustituyen- 
S,Í 0 
d Dn U 
o (17) en (16) se obtiene: U l )= Rb (18) 
uN p 0 y % 105) 
le + Se * ple). de 
s1% 
En esta última fórmula: Us, =lim, .,U, (e, )= (R, + R¿): Ys (19) 
» 


Así pues, las reobases para pulsos rectangulares de voltaje y las reobases para pulsos 
rectangulares de corriente se conectan entre sí a través de la resistencia total R 6 ik 


vista por el generador de pulsos en condiciones de corriente continua. 

El factor reobase de corriente en la fórmula (19) contribuye a incrementar la reobase de 
voltaje con el aumento en el radio del electrodo, en el espesor del casquete de tejidos no 
excitables interpuesto entre el electrodo y el miocardio y con el cambio en la modalidad 
de estimulación de monopolar a bipolar, puesto que la reobase es proporcional a 
di + y? (1, +e)?. No obstante en el otro factor presente en la mencionada fórmula, apa- 


rece R p ue se comporta en función del radio Y, aproximadamente como en . Puesto que 
» 

un factor aumenta y el otro disminuye al aumentar el radio del electrodo, cabe esperar que 

la reobase de voltaje presente un mínimo para un cierto radio. Este comportamiento de la 

reobase de voltaje, sobre el cual Werner Irnich llamó la atención por primera vez, se ha 

podido confirmar en trabajos experimentales [12]. 
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Si los electrodos son impolarizables, K 6 ) y Lle ) se anulan al igual que pl: ). En ese 
caso las curvas amplitud umbral versus duración se comportan igual para pulsos de co- 
rriente y para pulsos de voltaje. Pero si los electrodos se polarizan, de la fórmula (18) se 
desprende que cuando la duración del pulso aumenta, las amplitudes umbral de los pulsos 
rectangulares de voltaje se aproximan a su reobase más rápido que las amplitudes umbral 
de los pulsos rectangulares de corriente se aproximan a su propia reobase. 

Por definición, las constantes de tiempo para los pulsos de voltaje y de corriente vienen 


lim, _p(t, -U,[t,)) lim, o(£, +2, (e, )) 
Y yo E 
O Lio 


dados por los límites [19], [20] e 


Entonces de (5) y (18) resulta, teniendo en cuenta que pl0) Ea Rs: 
R 
Pp 


Loa R 
su = E 2 Ls (0) 
1+p(0) R,+R, 
La Figura 5 muestra algunas correlaciones entre el radio del electrodo y las constantes de 
tiempo del sistema electrodo-miocardio. 


ms 


” 


Figura 5. Resultados experimentales sobre la variación de las constantes de tiempo /;s , y 


É ¿ , con el radio 7, del electrodo de marcapaso. Tomada de la referencia [20]. 


Los resultados experimentales confirman las predicciones que se desprenden de la fórmu- 
la 20) para los electrodos convexos convencionales [20]. 

En igualdad de las demás condiciones, si los electrodos se polarizan, la constante de 
tiempo para los pulsos umbral de voltaje de forma rectangular es siempre menor que la 
constante de tiempo para los pulsos umbral de corriente de forma rectangular, y la dife- 


rencia entre ellas aumenta cuando aumenta el valor de la resistencia de la interfase R, en 


condiciones de corriente continua. 
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Además, para un intervalo de radios de electrodo tal que R » es un orden de magnitud 


mayor que R, ,ambas constantes de tiempo aumentan con el radio en forma análoga. 


(6) Características de una cabeza de electrodo de marcapaso que influyen sobre el 
umbral global de estimulación y sobre la calidad de la detección de la actividad eléc- 
trica del miocardio. 


Con el propósito de visualizar el umbral de excitabilidad global de un sistema electrodo- 
miocardio, la duración del pulso se grafica en abscisas y la amplitud (de voltaje o de co- 
rriente) se grafica en ordenadas en un cuadrante del plano. 

Si la curva amplitud umbral versus duración de los pulsos correspondiente a un primer 
sistema electrodo-miocardio se ubica por encima de la curva umbral correspondiente a un 
segundo sistema electrodo-miocardio, se puede decir que la excitabilidad global de este 
último es mayor que la excitabilidad global del primero. 


Para que esto ocurra es necesario y suficiente que se cumplan dos condiciones [22]: 


(1) La reobase R del primer sistema (sistema 1) debe ser mayor que la reobase del se- 

gundo sistema (sistema 2): R >R, 

(2) El producto de la reobase por la constante de tiempo? (o la cronaxia) del primer siste- 

ma debe ser superior a ese producto para el segundo sistema: —R,-f5,>R,-*ts, 

Si R,>R, pero R¡:tg,<R,:fg, hay una duración del pulso estimulador 
R) tg R,:t 


= 52 tal que si tp >1 
R,=R; 


el segundo sistema electrodo-miocardio es 


lo. P,x 


más excitable que el primero y si f, <f, el primer sistema electrodo-miocardio es más 


excitable que el segundo. 


Un objetivo de diseño es aumentar la excitabilidad global del sistema electrodo miocardio 
actuando sobre las características del electrodo. 

Las características de una cabeza de electrodo que se pueden modificar son las siguientes: 
- La macro estructura (tamaño y forma). 

- Los materiales. 

- La textura (micro-estructura) de la superficie. 


5 En el caso de umbrales de corriente, de la fórmula (8) resulta que el producto de la reobase por la 
constante de tiempo es igual a la carga umbral límite. Esto se utiliza en la sección 7 al comparar la 
excitabilidad global del sistema electrodo-miocardio de electrodo cóncavo-clotoidal con la corres- 
pondiente excitabilidad del sistema electrodo-miocardio de electrodo convencional. 
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-La capacidad de liberación lenta de fármacos que limiten la reacción inflamatoria local y 
minimicen el incremento fisiológico del umbral. 

En igualdad de las demás condiciones (fijación del electrodo por medios mecánicos inte- 
grados en el catéter, complacencia del catéter, estado del miocardio excitable en la zona 
adyacente al electrodo), la macro y micro estructura de la cabeza, así como la selección 
del material y la capacidad de liberar esteroides, deben permitir: 

-Umbrales globales, tanto agudos como crónicos, tan bajos como sea posible. 

-Un incremento fisiológico mínimo del umbral global que aparece asociado con la res- 
puesta inflamatoria del tejido miocárdico a la implantación del electrodo. 


Además de requerimientos relativos al umbral, es necesario contemplar requerimientos 
relativos a la detección de la actividad eléctrica del miocardio. La macro y micro estructu- 
ra del electrodo, así como la selección del material, deben conducir a las máximas ampli- 
tudes y velocidades de ascenso (slew rates) que sea posible para las señales eléctricas que 
ingresan a la unidad de estimulación y detección, mientras que las interferencias combi- 
nadas, debidas a interferencias de señales externas al sistema electrodo-miocardio deben 
ser mínimas [5], [6]. 


Como se mencionó en las secciones 4 y 5, dentro de ciertos límites, una disminución del 
área de la cabeza conduce a una disminución del umbral global: en igualdad de las demás 
condiciones la excitabilidad global del sistema electrodo miocardio se incrementa cuando 
el electrodo se hace más pequeño. Además, una disminución en el área de la cabeza del 
electrodo, dentro de ciertos límites, conduce a incrementar la duración de la batería del 
marcapaso, porque para una duración de pulso dada, al disminuir la amplitud umbral del 
pulso de corriente disminuye la carga eléctrica drenada en cada pulso. 


Un aumento del área de la cabeza mejora la calidad de la detección de la actividad eléctri- 
ca del corazón y disminuye los efectos de la polarización de la interfase electrodo-tejidos 
(polarización del electrodo). Para incrementar el área de la interfase sin aumentar el ta- 
maño del electrodo, se introdujeron los electrodos porosos [6], [27], [30]. 

Una disminución en la polarización de los electrodos es ventajosa cuando se aplican pul- 
sos rectangulares de voltaje, porque la fase descendente del correspondiente pulso de 
corriente conjugado (Figura 4 (b)) se hace menos marcada (el pulso de corriente se apro- 
xima más a un pulso rectangular), lo cual (como se dijo al final de la sección 3) permite 
alcanzar el umbral con amplitudes de pulso de voltaje menores. 


Las exigencias de desempeño en la estimulación y en la detección de la actividad eléctri- 
ca del miocardio, que están en buena medida contrapuestas, deben ponderarse en el dise- 
ño de los electrodos. A continuación se presenta un resumen de diseño conceptual de un 
tipo de electrodo cóncavo que intentaba contemplar esas exigencias en el marco de la 
tecnología y los conocimientos disponibles sobre el umbral en la primera mitad de la 
década de 1980. 
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(7) Un ejemplo de electrodo cóncavo: el electrodo clotoide. 


La idea de utilizar electrodos cóncavos apareció a comienzos de la década de 1980 del 
siglo pasado, en forma casi simultánea, en Francia y en Uruguay, cuando se había hecho 
evidente que hay una relación directa, no solamente entre la forma de la cabeza del elec- 
trodo y la estabilidad de su inserción en los tejidos, sino también entre dicha forma y las 
variables del umbral [6, 23, 31-31]. La Figura 6 muestra un esquema de la versión france- 
sa, el electrodo parabólico ASTER-1263. 
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Figura 6. Electrodo parabólico cóncavo-cóncavo ASTER-1263. 


En 1985 Amiel, Besson y Lehmann presentaron un trabajo donde resumen los fundamen- 
tos teóricos del diseño (lograr una convergencia de líneas de corriente sobre el miocardio 
adyacente), los cálculos y los resultados de pruebas clínicas de comparación de un elec- 
trodo convencional considerado de alta performance para la época, con el electrodo para- 
bólico que proponían [33]. Si bien el nuevo electrodo presentaba umbrales más bajos en 
agudo, las diferencias entre los umbrales de los electrodos que se compararon desapare- 
cían en estado crónico. Estos resultados condujeron a una pérdida de interés en el electro- 
do parabólico. 

La idea de emplear electrodos cóncavos en marcapasos cardíacos se le ocurrió al autor* en 
1983 mientras estudiaba los efectos mecánicos del flujo de agua en medios porosos de 
matriz granular [34]. Debido a una iniciativa del Dr. Orestes Fiandra”, desde 1979 el 
autor había estado investigando, desde el punto de vista físico-matemático, los problemas 
que plantea el umbral de estimulación eléctrica del miocardio. 


El Profesor Fiandra había inventado, construido e implantado un electrodo de auto- 


6 Desconociendo lo que ocurría en Francia. Se enteró del electrodo parabólico cóncavo-cóncavo 
recién en 1985, a través del Dr. Orestes Fiandra que había regresado de un congreso con un libro de 
actas donde Amiel, Besson y Lehmann presentaron los resultados de sus investigaciones. 


7 Orestes Fiandra fue un uruguayo eminente y exitoso: médico cardiólogo, catedrático de medicina, 
inventor y empresario, creador del Centro de Construcción de Cardioestimuladores (CCC). 
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fijación cuyo área de material conductor en contacto con el miocardio era mucho mayor 
que el área correspondiente de los electrodos convencionales que se utilizaban en aquella 
época, pero que inexplicablemente presentaba umbrales más bajos, dando origen a siste- 
mas electrodo-miocardio más excitables [32], [35]. 

Pensó que para comprender estos resultados al parecer paradójicos era necesario un enfo- 
que matemático del umbral. Por este motivo invitó al autor a participar en la investigación 
de éste y otros problemas que se planteaban a propósito de los electrodos de marcapaso. 


Hacia fines de 1982 los resultados de las investigaciones sugerían que un electrodo de 
marcapaso debería ser capaz de: 

(a) Generar una región crítica, de un tamaño y forma adecuados, donde las densidades de 
corriente fueran más intensas y relativamente paralelas, para abatir los umbrales de esti- 
mulación del miocardio. 

(b) Saltar por encima del casquete de tejidos excitables para contrarrestar el incremento 
fisiológico del umbral. 

(c) Disminuir lo más posible la densidad de corriente eléctrica en la superficie del elec- 
trodo, para una corriente total prefijada, minimizando la polarización de la doble capa 
electroquímica en la interfase entre el electrodo y los tejidos. Además, el electrodo debe- 
ría ser capaz de detectar con fidelidad la actividad eléctrica del miocardio. 


Lo bastante cerca de cualquier superficie regular de un electrodo las líneas de corriente 
son paralelas entre sí. Lo bastante cerca de una zona convexa con curvatura más elevada 
la densidad de corriente es más intensa y los umbrales en estado agudo pueden ser meno- 
res, como ocurría con el electrodo de auto-fijación mencionado. Pero el efecto de la cur- 
vatura en los electrodos convexos se pierde cuando se produce la reacción inflamatoria al 
cuerpo extraño, una vez implantado el electrodo en forma crónica [16], [17]. 


La idea de utilizar un electrodo cóncavo se planteó con el fin de subsanar estas desventa- 
jas de los electrodos convexos. La Figura 7 representa gráficamente lo esencial de la idea. 


Figura 7. La curva cóncavo-convexa en el corte indica la superficie del electrodo. La región 
crítica del electrodo se ha esbozado con una mancha en tonos de gris, junto con algunas líneas 
de corriente y trazas de equipotenciales. Se puede ver el efecto de las curvaturas en los extre- 
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mos convexos de la interfase y el relativo paralelismo y concentración de líneas de corriente en 
la región crítica situada sobre el eje del electrodo. El esquema es cualitativo. 

A partir de esta idea básica, el autor exploró en forma preliminar la posibilidad de diseñar 
una clase de electrodos cóncavos con simetría de revolución (meridiano-plana) minimi- 
zando los efectos negativos de las inevitables zonas convexas situadas en los bordes de la 
concavidad y maximizando los efectos positivos asociados con la región crítica. 

Este trabajo condujo a concebir una clase de electrodos cóncavos cuya superficie interse- 
ca un semi-plano por el eje de simetría del electrodo en una curva que se construye a 
partir de un tipo de curva en espiral, bien conocida en ingeniería vial como clotoide o 
espiral de Cornú y en óptica física en relación con las integrales de Fresnel de la teoría de 
la difracción (las funciones integrales de Fresnel se emplearon en el diseño). 


Luego se llevaron a cabo investigaciones experimentales en una cuba electrolítica con 
siete modelos de electrodos a escala, dos convexos y cinco cóncavos (de los cuales dos 
eran clotoides), entre fines de 1983 y comienzos de 1985. Los resultados obtenidos sugie- 
ren que, efectivamente, para un electrodo cóncavo se forma una región a cierta distancia 
de su superficie, donde la densidad de corriente es más elevada en comparación con los 
electrodos convexos del mismo tamaño y material, y las líneas de corriente son relativa- 
mente paralelas (región crítica del electrodo). 

Además, esta región crítica se encuentra más alejada de la superficie del electrodo y es 
más grande para los electrodos clotoides en comparación con los otros electrodos cónca- 
vos investigados, del mismo tamaño y material [36]-[39]. 

El Dr. Fiandra y el autor seleccionaron un electrodo clotoide para implantar en un modelo 
animal junto con un electrodo convexo convencional del mismo tamaño (2 mm de diáme- 
tro) y material (platino). La Figura 8 (a) muestra la plantilla que preparó el autor. 

A partir de ella se hizo (en Barcelona) una herramienta de corte especial que el Dr. Fian- 
dra utilizó en su torno de joyero para producir la concavidad en un cilindro de platino. 

La Figura 8 (b) muestra una fotografía de un electrodo clotoide ya terminado (incluyendo 
el proceso que lo transforma en un electrodo poroso). 


(a) (b) 


Figura $. A la izquierda se ve la plantilla con la curva seleccionada. A la derecha se ve un 
electrodo ya terminado. La regla da una idea de las dimensiones de la concavidad y del diáme- 
tro del electrodo. 
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Entre 1985 y 1986 se efectuaron experimentos de tres semanas de duración con modelos 
animales utilizando un electrodo plano y uno clotoide en el mismo corazón. 

Entre 1987 y 1993 se patentó la invención en Uruguay [36] y se efectuaron pruebas clíni- 
cas en España. 

A continuación, a vía de ejemplo, se presentan valores representativos de las reobases de 
corriente, las constantes de tiempo correspondientes y las cargas umbral límite obtenidas 
en dos modelos animales (D-4 y D-5) [37]. 

El electrodo clotoide se representa mediante (a) y el convencional mediante (f). 

Los datos corresponden al día en que se implantaron (Día 0), al día siguiente (Día 1), al 
cuarto (Día 4) y al octavo día (Día 8) a partir de la implantación de los dos electrodos en 
el mismo corazón. 


Modelo D-4 Día 0 Día 1 Día 4 Día 8 
Reobase (0) 0.20 0.71 0.71 0.72 
(mA) B 0,49 0.68 1.12 1.13 
Constante de 
tiempo (0) 0,54 0.36 0.21 0.38 
(ms) B 0.35 0.15 0,29 0.30 
Carga umbral 
límite (0) 0.11 0.25 0,22 0.27 
B 0.17 0.10 0.32 0.34 
Modelo D-5 Día 0 Día 1 Día 4 Día 8 
Reobase (0) 0.64 1.22 0.84 0.81 
(mA) B 0.39 1.05 1.55 2.15 
Constante de 
tiempo (0) 0.27 0.21 0.17 0,49 
(ms) B 0.21 0.24 0.15 0.31 
Carga umbral 
límite (0) 0.17 0.26 0.14 0.40 
B 0.08 0.25 0.23 0.67 


Para comparar el umbral de corriente del electrodo clotoide (a) con el umbral del electro- 
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do convencional ($) se comparan reobases y cargas umbral límite como se describió en la 
sección 6: una reobase menor junto con una carga umbral límite menor corresponden a un 
umbral menor (o sea, una excitabilidad global mayor). 


En general, en el momento del implante la curva amplitud umbral versus duración de 
pulso para el electrodo clotoide se dispone ligeramente por encima de la curva correspon- 
diente al electrodo plano, pero luego de transcurridos varios días la situación se invierte 
porque las reobases y las cargas umbral límite del electrodo convencional aumentan en 
forma significativa (incremento fisiológico del umbral) mientras que las reobases y las 
cargas umbral límite de los electrodos clotoides aumentan mucho menos, mantienen sus 
valores en agudo o incluso disminuyen. 


Las constantes de tiempo de ambos tipos de electrodo pueden disminuir primero respecto 
de los valores que presentan en el momento del implante, para luego aumentar y estabili- 
zarse. 


Por último, en general el electrodo clotoide mostró una mejor performance en la detec- 
ción de la actividad eléctrica del miocardio respecto del electrodo convencional. En las 
referencias [21], [26] se puede hallar una investigación sobre las propiedades de detección 
de señales por electrodos cóncavos. 


(8) Conclusiones 


-Para estudiar los problemas que se plantean en relación con los umbrales de estimulación 
eléctrica del miocardio parece conveniente introducir una distinción entre excitabilidad 
global y excitabilidad local en el marco de una descripción funcional del sistema electro- 
do-miocardio como la propuesta en la sección 2. 


-Para estimulación mediante pulsos rectangulares de corriente controlada, las fórmulas 
obtenidas en la sección 4 predicen el aumento de las reobases de corriente y las constante 
de tiempo para pulsos de corriente, con el incremento del radio del electrodo, el incre- 
mento del espesor del casquete de tejidos no excitables interpuesto entre el electrodo y el 
miocardio y con el cambio en la modalidad de estimulación de monopolar a bipolar. 


-Para estimulación mediante pulsos rectangulares de voltaje controlado, las fórmulas 
obtenidas en la sección 5 relacionan las reobases de voltaje y las constante de tiempo para 
pulsos de voltaje con la resistencia variable de los tejidos, la resistencia de la capa doble 
electroquímica de la interfase a la corriente continua, y con las reobases de corriente y las 
constantes de tiempo para pulsos de corriente, respectivamente. Predicen las modificacio- 
nes en el umbral de voltaje con las variaciones en las dimensiones, forma y polarización 
del electrodo, así como con las modificaciones en la resistencia de los tejidos asociadas a 
la formación del casquete de tejidos no excitables. Si bien la evidencia experimental co- 
nocida por el autor confirma todas las predicciones, en general esa confirmación es cuali- 
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tativa. Esto es debido a que la evidencia proviene de experimentos que no se diseñaron 
especificamente para poner a prueba las fórmulas presentadas y discutidas en este trabajo. 


-En el fenómeno de histéresis que presenta el umbral de estimulación del miocardio, 
además de una componente puramente electrofisiológica, relacionada con el denominado 
efecto Wedensky, hay una componente mecánica relacionada con el comportamiento del 
miocardio como un sistema mecánico sobre-amortiguado y sometido a las fuerzas perió- 
dicas producidas por la contracción del corazón. El aporte de la componente mecánica a 
la histéresis del umbral es un tema para investigar. 


-S1 se admite que la interfase electrodo-tejidos se puede describir mediante una impedan- 
cia de Drake, entonces una forma de disminuir la polarización de un electrodo de marca- 
paso, contribuyendo a disminuir el umbral de voltaje en el sistema electrodo-miocardio, 


es disminuir los parámetros M y K, que aparecen en la fórmula (10 c). Puesto que 


K, ¿es inversamente proporcional al área de la interfase, conviene aumentar esta última 


sin aumentar el área aparente (el radio del electrodo). 
Como m<3-—d, (fórmula (10 b)), se puede disminuir el parámetro 1M aumentando lo 


más que sea posible la dimensión fractal d s de la interfase desde 2 que corresponde a 


una superficie idealmente lisa, hacia 3 que correspondería a un fractal que llena el espacio 
(lo cual, por supuesto, no es alcanzable en la práctica). En general los procedimientos que 
aumentan la porosidad de un electrodo incrementan simultáneamente su dimensión frac- 
tal. 


-A diferencia de lo ocurrido con el electrodo cóncavo-parabólico francés, al final del 
período de tiempo abarcado por los experimentos la excitabilidad del sistema electrodo 
miocardio formado por el electrodo cóncavo-clotoide resultó mayor que la del electrodo 
convencional utilizado para la comparación. Una de las razones de esta diferencia de 
comportamiento entre los dos electrodos cóncavos, además de la diferencia en la forma de 
la concavidad, es la diferencia en las dimensiones de los electrodos: el francés es bastante 
más pequeño [1] lo cual desvanece el efecto de la concavidad cuando se forma el casquete 
de tejidos no excitables. 


-Los avances en la tecnología de generación de pulsos estimuladores y en la tecnología de 
detección de señales en los marcapasos han sido mucho más rápidos que el avance en el 
diseño de los electrodos. Un motivo de este retraso parece ser una comprensión insufi- 
ciente de lo que ocurre en el miocardio cuando sus membranas excitables se polarizan 
debido al campo de densidades de corriente generado por un electrodo externo. 
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